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ABSTRACT 

Solutions are derived for a class of boundary value 
problems for the time -dependent temperature distribution in 
a two layer, composite slab with contact resistance at the 
interface and contact or convective resistance on the outer 
boundaries. The results of a parametric computer study 
using these solutions are presented. This study includes a 
set of dimensionless correlations of an arbitrarily defined 
time to approach steady state and a discussion of some tran- 
sient thermal phenomena which are characteristic of systems 
of this type. 

The applicability of these solutions for the predic- 
tion of the transient behavior of real composite systems is 
evaluated by comparison with experimental results. Experi- 
mental data are presented for the transient thermal response 
of several composite metal systems when subjected to thermal 
transients which closely approximated the boundary condi- 
tions of the theoretical solutions. 

Comparisons with the experimental data indicate that 
the theoretical solutions can be used to predict the tran- 
sient thermal response of systems to which they are applica- 
ble to an accuracy sufficient for most engineering purposes. 
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I . INTRODUCTION 


Objectives and Scope 

It is the primary objective of this writing to pre- 
sent the results of theoretical and experimental investi- 
gations made by the writer on one -dimensional composite 
systems with contact resistances when subjected to thermal 
transients. These results are believed to be a significant 
contribution to the fields of heat transfer in composite 
media and thermal contact resistance. In the remainder of 
this section some of the basic definitions and ideas asso- 
ciated with the concept of thermal contact resistance are 
presented to form a background and illustrate the impor- 
tance of the problem of contact resistance. 

In section II the actual mechanism of contact heat 
transfer is examined and discussed with a view toward under 
standing what variables influence the phenomenon of contact 
heat transfer. 

In section III the existing literature is examined 
with a discussion of some of the major works in the field 
and a categorization of the literature by the type of infor 
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mat ion to be found in the various references. 

The derivations of the analytical solutions obtained 
by the writer for the time -dependent temperature distribu- 
tions in composite solids with contact resistances are giv- 
en in section IV. The results of a computer study made 
with these solutions are also presented and discussed.' 

Section V contains the descriptions of the equipment 
and procedures used in obtaining the experimental data. 

The experimental results are presented in section VI. Con- 
clusions drawn from a comparison of the experimental and 
theoretical results are presented in section VII. Recom- 
mendations for future investigations are also offered in 
section VII. 

Background 

The particular area of interest in heat transfer 
which is concerned with the transfer of heat across sur- 
faces in contact is relatively new. Truly active interest 
and investigation in this area did not begin until approx- 
imately eighteen years ago. Contact resistance studies are 
basically concerned with the effects on heat transfer rates 
of the presence of a joint or interface surface of contact 
between two solid bodies. Although the majority of the back 
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ground material in the field is concerned with steady state 
contact phenomena, it is presented here because of its ob- 
vious importance to the study of non-steady contact resist- 
ance. In order to introduce the subject it is logical to 
discuss how the surfaces of solid bodies make contact with 
each other. 

A surface on a solid body which has been formed by a 
machining process is usually thought of as smooth. There 
are degrees of "smoothness" for such a surface, depending 
on the nature of the machining process. To some extent the 
differences can be qualitatively distinguished by the human 
senses of sight and touch. However, even though a surface 
may appear to these senses to be quite smooth, it is known 
that on a microscopic level the surface is rough. That is, 
under sufficient magnification, a surface is not smooth but 
would appear as a series of irregular hills and valleys, 
much like an aerial view of a mountainous terrain. 

All machined surfaces possess some degree of rough- 
ness. Consequently, when two such plane^ surfaces are 
pressed together they can be in actual contact only at dis- 


Although most of general remarks apply to all ma- 
chined surfaces, the present work is concerned only with 
plane, i.e., nominally flat, surfaces. 
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Crete points. Figure 1 illustrates how the interface be- 
tween two solids might appear. The two bodies touch only 
where peaks touch peaks or peaks touch valleys. These con- 
tact points are variously called "contacts," "contact 
spots," and "a-spots," The latter designation was picked 
up by the early writers from the pioneer work in electrical 
contacts by Holm [ll8] . Thus the term a-spot is the most 

prevalent among writers in the field and it will be used 
in the present work. The total area of actual contact be- 
tween two solids is the sum of the areas of the individual 
a-spots. This total area may be less than 1% of the appar- 
ent contact area, and is above 20% only for carefully pre- 
pared surfaces [41, 43, 78] . 

The fact that the interface contact between two solids 
occurs only at discrete points gives rise to the heat trans- 
fer phenomenon known as "contact resistance." Such an in- 
terface, generally referred to as a "contact" or "joint," 
causes a constriction of the heat flux lines (Figure 1) . 

Thus an additional resistance to heat flow is produced be- 
cause the heat flow, which is parallel at some distance 


2 

Numbers in brackets refer to the corresponding ref- 
erences listed in the bibliography. 
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-Microscopic View of Two Surfaces in Contact 
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from the interface, must bend to pass through the a-spots. 
In treating this phenomenon it is a matter of choice as to 
whether one uses resistance or its reciprocal, the conduc- 
tance. Most of the authors in this field have chosen to 
deal with conductance. Therefore, the conductance will be 
used here (however, the term "resistance" will be used in 
certain instances because of its intuitive appeal) . 

Before examining the details of the heat transfer 
problem of a contact interface it is desirable to give the 
accepted definition of contact conductance and to comment 
on it. Contact conductance is defined (analogous to a con- 
vection heat transfer film coefficient) as follows: 



Ac UDc 


(Btu/hr-ft 2 -°P) (1-1) 


In which Q c = heat transfer rate across the interface, A c 
is the apparent contact area (cross sectional area of sol- 
id) , and (AT) C = the apparent temperature drop at the con- 
tact. This definition introduces the fiction of an "appar- 
ent" temperature drop at the interface. Obviously there is 
no real discontinuity of the temperature distribution 
through the solid contacts. There is a continuous distri- 
bution of temperature extending through the contact from 
both solids. However, such a temperature distribution would 
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be extremely difficult to describe. From Figure 1 it is 
clear that the distribution would be three-dimensional even 
if the distribution in main part of the bodies was one-di- 
mensional. For this reason the concept of a contact tem- 
perature drop is a convenience. This temperature drop is 
defined, as shown in Figure 2, as the difference in the 
temperature obtained by extrapolating the temperature pro- 
files in the two regions to the interface boundary. The 
description of this definition as fictional is not intended 
to imply disapproval. On the contrary, it is a practical, 
working definition. The only caution required here is that 
if one obtains the temperature profiles experimentally, care 
must be taken to insure that the measurements are made at a 
sufficient distance from the contact to insure that the iso- 
therms are essentially parallel planes. 

Thus the effect of a contact interface can be thought 
of either as a discontinuity in the temperature distribution 
on a macroscopic level, or as a constriction of the heat 
flow lines on a microscopic level. In either case the im- 
portant point is that the interface causes an additional re- 
sistance to heat transfer. The prediction of this resist- 
ance for a given set of conditions, or more commonly, its 
reciprocal, conductance, is the primary goal of all the the- 


DISTANCE 


Pig. 2. — Definition of Contact Temperature Drop 
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oretical and experimental work in this field. 

Applications 

The study of contact conductance is an important area 
in the field of heat transfer. Practically all types of 
equipment and hardware in which the transfer of heat is of 
importance contain composite structures, or joints, between 
members or pieces of equipment, etc. Each such joint that 
does not provide intimate contact, i.e., any joint that is 
not welded, brazed, etc., produces a contact resistance. 

The importance of the knowledge of contact phenomena is par- 
ticularly significant in situations in which the conductance 
may vary with time or over the range of operating conditions 
of the equipment. This is especially true if the contact is 
part of a critical system of heat removal or thermal isola- 
tion. A few such applications are: equipment for heat re- 
moval and environmental control for electronic components, 
high temperature heat -treatment equipment, nuclear reactors, 
high speed aircraft, satellites, and launch and reentry ve- 
hicles. The National Aeronautics and Space Administration 
is currently very interested in contact conductance - as ev- 
idenced by their sponsoring of several research projects in 
this field, including the present work. 
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Although much work has been done, both theoretically 
and experimentally, almost all of it has been done for 
steady-state conditions. These works have produced only a 
reasonable understanding of steady-state contact conduct- 
ance. As will be discussed in section III, there is a need 
for more work to be done in the study of steady -state con- 
ductance. On the other hand, very little has been done for 
non- steady conditions, and it is for this reason that the 
present work was undertaken. 


II. the mechanism of contact heat transfer 


The definition of contact conductance was given in 
Section I in equation (1-1). It was demonstrated by means 
of Figure 2 how one could calculate the conductance, he, 
from experimental data. However, in order to be able to 
predict values of h c for a given situation it is necessary 
to know what the parameters are that affect it. In other 
words, a knowledge of the physical mechanism is required. 

It is the purpose of this section to present a description 
of the physical mechanism of heat transfer across a contact. 
Such a presentation must include a description of the con- 
ditions under which certain variables exert a considerable 
influence and when their influence is negligible. Thus 
there may be some overlap between the conclusions stated 
here and the discussions in the next section. However, the 
small amount of repetition is felt to be justified for the 
sake of presenting an overall picture here. 

If there is a flow of heat between two solid bodies in 
contact the mechanism by which the heat flux traverses the 
contact plane is, in general, very complex. As illustrated 
in Figure 1, heat may be transferred across an interface 
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in the following ways. 

1) solid-to-solid conduction heat transfer through 

the contact points. 

2) conduction heat transfer through the interstitial 

gas (if a gas is present). 

3) solid-to-solid radiation heat transfer between 

the portions of the surface that do not touch, 

i.e., through the interstitial volume. 

4) convection heat transfer through the interstitial 

gas. 

Although the above heat transfer modes are somewhat 
interdependent, for the purposes of analysis they will be 
assumed to be independent. Such a simplification of con- 
ditions always creates some doubt about conclusions drawn 
on this basis. However, as will be shown, the resulting 
conclusions indicate that the differences in the amounts of 
heat transferred by the different modes are orders of magni- 
tude apart. Thus, the conclusions are believed to be justi- 
fied since it would be difficult to imagine that the sim- 
plification of the model would produce distortions of this 
size. With the assumption of independence the total con- 
ductance, h, of the contact interface can be assumed to be 
the sum of the four individual conductances: 
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+ + fir f (II-l) 

where the subscripts refer to solid, fluid, radiative, and 
convective conductances, respectively. 

Figure 3 is an idealized representation of a contact 
interface which will be used to analyze the importance of 
the various contact heat transfer modes. 

Interstitial Convection 

It will be assumed that the contact points are widely 
spaced, i.e., the size of the contact spots (their radii) 
is small compared to the distance between them, see Figure 
3. It is also assumed that the gap thickness, 5 { , is small 
compared to the distance between contact spots. With these 
assumptions the convection analysis can be based on two 
large parallel flat plates separated by a distance Sf (the 
effective fluid gap thickness) . Three different cases of 
orientation could arise: 1) gap horizontal with the upper 

plate at a higher temperature; 2) gap horizontal with the 
lower plate at a higher temperature; and 3) gap vertical. 
For the first case no convection will occur [245, p. 272], 

For the second and third cases, experimental work has shown 
that no convection occurs unless the Grashof number (based 


Heat Channel 
With Radius *' 
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Fig. 3. — Idealized Contact Geometry 
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on gap thickness) is greater than 1700 and 2000, respec- 
tively [245, p. 272 —3] . In the type of engineering surfaces 
normally encountered the Grashof number (based on and 
properties at the mean interface temperature) ranges from 
10 ^ to 10 Therefore the convective heat transfer for 
any realistic contact problem is safely ignored. 

Gap Radiation 

With the same assumptions made above, the radiative 
heat transfer analysis can also be based on large parallel 
flat plates. For such a configuration the radiative heat 
transfer rate per unit area is given by Jakob [246, p. 3] as 

$'-•* 1 £ x . <r ( t , «- ( n - 2 ) 

where r is the Stefan -Boltzmann constant, and T 2 are the 
respective surface absolute temperatures, and 6- ^ and £ 2 
are the emissivities of the surfaces. If the assumption of 
gray surfaces is made the interchange factor is given by 

= Ye, + Vt x -7 (II " 3) 

Using (II —2 ) and the conductance definition the fol- 
lowing relation for radiative conductance results 


r (r,-r 2 ) 


( II -4) 
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or, 

^.= £ /x <s- (T’+r/Tt+^rS+T;*) di-si 

If the assumption is made that the temperature differ- 
ence (T^ - T 2 ) is small compared to the mean temperature 

(T « T 1 + T 2) , then (II-5) becomes 

m 

2 

For most thermal contacts the above assumptions are 
reasonable. It can be seen from (II -6) that the contact 
conductance would depend strongly on the mean interface tem- 
perature if radiation were a considerable portion of the to- 
tal heat transfer. Experimental work has shown that this 
dependence does not occur (24, 27, 47, 140] . There is some de- 
pendence on T m , but the extent to which it appears indicates 
that it is primarily the result of the dependence of the 
solid and fluid conduction modes. Using (II-6) and exper- 
imental results for the total conductance, Fenech and 
Rohsenow [84] have found that for mean contact temperatures 
below 1100°F the radiative heat transfer amounts to less 
than 1 % of the total. Similarly, Clausing [58] has stated the 
radiation accounted for less than 2 % in the worst conditions 
of his experimental work. Therefore except in the presence 
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of very high temperatures or very poor conductors the con- 
tribution of gap radiation to contact heat transfer is neg- 
ligible. In the present experiments the gap radiation con- 
ductance coefficient was estimated from (II-6) to be about 
0.2 Btu/hr.-ft?-°F. Comparing this with lowest value of 
the total conductance coefficient measured, 29 Btu/hr.-ft.- 
°F, shows that it is less than 1 %. 

Thus, in the remainder of this work, the terms hg and 
h r in equation (II-l) are assumed to be negligibly small. 
The total conductance will be written unsubscripted and the 
solid and fluid components will be subscripted with ”s" and 
"f" respectively, i.e., 

h = total contact conductance = h + h, 

s f 

Gas Conduction 

Utilizing the same assumptions about the gap size and 
spacing as were used above, the interstitial fluid conduct- 
ance can be written as the reciprocal of the fluid resist- 
ance , 

where is the "ef fective" gas thickness and 5^- is the 
fluid thermal conductivity. The prevailing theories and 
practices regarding the evaluation of are presented in 
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the next section. It suffices here to say the is small 
and it is related to the surface roughness. By means of 
equation (II -6) it can be seen that because $£ is small, h^ 
can be quite large even though k is typically small for 
most gases and some liquids. In some cases even with air in 
the interstices the fluid conductance may account for a ma- 
jor portion of the heat transfer. For example, Barzelay, et 

al [24, 25] report total conductances as low as 250 Btu/hr.- 
2 o 

^ f° r a stainless steel contact with a surface rough- 
ness of 120 microinch at a mean interface temperature of 
200 F (for low contact pressures where solid conduction is 
small). For such a case k^is approximately .0181 Btu/hr. - 
ft.- F [245] and S f is approximately 3.0 x 10~ 4 ft. With 
the values equation (II -7) gives 

£ = — sr to a r*. 

■* 3* jo' 9 

This would indicate that for this case the fluid conductance 
accounts for approximately 25% of the total heat transfer. 
Cases have been reported where the fluid conductance ac- 
counts for over half the total [41] , and even as high as 98% 
0-37] . Equation (II -7) also brings out another important 
point: that one can accomplish a reduction in contact re- 

sistance by putting a fluid with high thermal conductivity in 


the interstices. 


For most gases with the pressures normally encountered 
the dependence of thermal conductivity on pressure is small. 
However at low pressures where the mean free path of the gas 
molecules is of the order of, or greater than, the gap 
thickness this no longer holds. Since thermal contacts in 
a vacuum are of considerable importance, especially in the 
present work, the influence of low pressures on the gas con- 
ductivity will be discussed here. 

Again referring to Figure 3, the analysis will be 
based on two large flat parallel plates. For this configu- 
ration the following equation for the heat transfer rate 
per unit area is given by Kennard [247 p. 3173 » 



cl, a 


L 





(II-8) 


In which a^ and are the accommodation coefficients 
hhe two surfaces, ^ , R and c are gas density, gas 

constant, and constant -volume specific heat, respectively. 

convenient to define an "effective" accommodation co- 
efficient as follows, 

a = — — 
a t +ct 2 -cL t ci 2 

Equation (II-8) can be put into the form in which it custom- 
arily appears by using the above definition and the follow- 



ing thermodynamic relations: 

P . 6 

f = er 1 Cv ■ *-/ 

Where is the ratio of specific heats. Thus (II-8) be- 
comes , _ - i 

with the definition of h^ this becomes 

** a (II - 9) 

The prime notation was introduced to account for accommoda- 
tion effects. The T* is the mean of T£ and T^ which are 
respectively, the temperatures corresponding to the mean 
speeds of the molecules leaving the respective surfaces. 

The p* is the pressure of a gas which has the density of the 
gas in the gap but which is maxwellian and at a temperature 
T* . The T ' can be calculated from the actual plate tempera- 
tures and the accommodation coefficients. It is sufficient 

for present purposes to realize that if a and a do not 

1 2 

differ greatly, T‘ will be close to the mean of T and T 

1 2 

(if = a 2 * T is exactly T^) . Thus for most commonly 

encountered contacts p' and T' may be taken as p and T . 

m 

With these statements and equation (II-9) it is readily seen 
that the free molecule fluid conductance is inversely pro- 
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portional to the square root of the mean temperature and 
directly proportional to the pressure. The significant 
point about (II -9) is that the conductance is independent 
of the gap thickness. It should be emphasized that this 
conclusion only holds for the conditions of the derivation 
(II-8), namely, free molecular motion. The usual crite- 
rion for free -molecular motion is that the Knudsen number, 
v? f , be greater that 10, where \ is the mean free mole- 
cular path. In most of the published contact conductance 
investigations the maximum value of is less than .002 

inches. Using the above criterion and this one finds 

that equation (II-9) is valid for mean free paths greater 
than .02 inches, or, in terms of pressure, for pressures 
less than 100 microns of H . The most common values of 
that are usually encountered are less than the above and 
thus the assumption of free molecular flow is valid for low- 
er mean free paths, i.e., for higher pressures. In the 
present experimental work the largest Sf was calculated to 
be less than .0002 inches. 

It has been shown how one could calculate the value 
of the fluid conductance for the conditions of normal pres- 
sures and very low pressures. The former requires the de- 
termination of the effective gap thickness, which is usually 
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not a simple task, whereas the latter requires a knowledge 
of the mean gap pressure and temperature, the gas proper- 
ties and R, and the accommodation coefficients. For a 
situation in between these extremes it would be very diffi- 
cult to evaluate gas conductivity analytically. Fortunately, 
however, the conditions of interest are usually either at- 
mospheric or high vacuum conditions. 

Solid Conduction 

For most cases of practical importance the solid con- 
duction mode accounts for the largest part of the total con- 
tact heat transfer. Exceptions have been noted in the pre- 
vious paragraphs. Because solid conductance is the predomi- 
nant mode, especially in a vacuum, it has received a great 
deal of attention. 

In the absence of oxide films or other surface contam- 
ination this mode of contact heat transfer can be treated as 
simple conduction. The additional resistance to the flow of 
heat is caused by the fact that the heat flux lines must con- 
verge to pass through the areas of solid-solid contact. This 
constriction of heat flow lines has led to the use of the 
term "constriction resistance" or "constriction conduct- 
ance". The analysis of the solid conduction contribution 
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may be carried out in several ways. In the following para- 
graphs four separate methods of calculating the solid con- 
ductance term are presented. These four represent the ma- 
jor contributions thus far published. Though there have 
been other works published, they are only adaptations or 
modifications of the ones presented below. 

There are some differences and some similarities be- 
tween the four approaches which should be emphasized here. 
For this purpose it is necessary to understand the two kinds 
of surface irregularities which can exist on a solid surface. 
All surfaces possess a certain roughness. A surface may al- 
so have waviness. Surface waviness is the macroscopic "non- 
flatness" of a surface as opposed to surface roughness which 
is microscopic. Figure 4 is a sketch which shows how these 
surface characteristics might appear. Obviously the heat 
flowing across a contact formed by bringing together two 
surfaces which have both roughness and waviness could suffer 
two types of constriction. The heat flow lines would have 
to converge first to the "macroscopic contact areas" which 
exist because of the waviness. Then, within these areas, 
they would have to converge again to pass through the micro- 
scopic contact spots. Thus the terms "macroscopic constric- 
tion resistance" and "microscopic constriction resistance" 


Roughness Portion 



Fig. 4. — Surface Roughness and Waviness 
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have come into use. The basic difference between the four 
a PP roa< ~hss to be discussed lies in the assumptions made re- 
garding which of these two constriction resistances is pre- 
dominant. The first three of the four approaches which are 
discussed below are similar in that they assume that the mi- 
croscopic resistance dominates the contact resistance. They 
differ however in that the first two start with an attempt 
to solve for the temperature distribution in the immediate 
vicinity of the contact, whereas the third uses a simpler 
a PP r °ach . Two of three do employ the concept of waviness 
in some way, whereas the other assumes that the microscopic 
contact spots are distributed uniformly over a non-wavy sur- 
^ a ^ e • The fourth method assumes that the macroscopic con- 
striction resistance dominates. All four methods are dis- 
cussed individually below. A brief discussion of the sig- 
nificance of their basic assumptions and the limits of their 
applicability is presented at the conclusion of this section. 

The Approach of Cetinkale and Fishenden: 

Cetinkale (Veziroglu) and Fishenden [55] made a theoret- 
ical analysis of thermal contact conductance using a model 
as depicted in Figure 3. In this model there are N such con- 
tact spots over the entire surface. The heat flow at a large 
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distance from the contact interface is assumed to be paral- 
lel and uniform throughout the solid. Figure 5 is a sketch, 
to a more realistic scale, of the idealized contact element 
showing approximately how the isotherms and heat flow lines 
should appear. Also shown is the "dividing flow line" which 
separates the heat flow through solid spot from that which 
flows through fluid gap. The isothermal surfaces in the vi- 
c inity of the solid spot were assumed to be ellipsoids of 
revolution and the flow lines are then hyperbolas. An ex- 
pression for an elemental resistance was written in terms of 
the two parameters which occur in the equations of the el- 
lipses and hyperbolas. The expression was integrated to 
yield an equation for the contact resistance in which only 
one parameter remains to be eliminated. The authors deduced 
the value of the remaining parameter from the results of re- 
laxation solutions for a series of cases which covered a 

representative range of quantities a/b and bk /£ k 

f w f m 

After some further manipulation the following expression for 
the solid contact conductance is obtained. 

(f)*~ 

c - /? 

f > 5 <*/ 

It can be seen that the resulting expression is a transcen- 
dental equation for h and thus must be solved by some tri- 
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al and error procedure. Of the parameters which appear in 

equation (11-10) , k (harmonic mean of k and k ) should 

m 12 

be known for a given set of materials. For cases in which 
radiation can be neglected and the fluid is either a gas at 
higher pressures or a liquid, the effective conductivity of 
the gap is simply the fluid thermal conductivity. Otherwise 
the quantity k^ / S w hich is h , can be calculated by 
combining the results of the previous sections of the pres- 
ent work. Cetinkale (Veziroglu) and Fishenden call the ra- 
tio a/b the "constriction number." They assume that when 
pressure is applied to the contact, the softer of the two 
contact metals flows plastically until the average a-spot 
pressure is equal to the Meyer hardness. Thus, since 





Atast - 

Ajw 6 X 


and since 



Apparent Contact Area 
Actual Contact Area 


(II-ll) 


pA = MA = total applied load 
c 


the constriction number can be calculated from 

% "YW (ii ' i2) 

For those cases where 8 f is needed separately, the authors 


suggest that 




(H-13) 
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Where e is a geometric factor depending on the shape of 
the surface irregularities and ^ and are the respec- 

tive arithmetic mean heights of the surface roughnesses. 

The authors noted that for ground surfaces, $ is almost 
constant. Finally, the heat channel radius, b, must be 
found from the relation 


4 - r-u+u-ar 


(11-14) 


Where \ ^ and A are the wavelengths of roughnesses of 
the two surfaces. The quantities ^ and V are constants 
which must be experimentally determined using heat transfer 
results. Experimental results obtained by these authors for 
ground surfaces were: 0.61, 0.0048, and - 5/3 for e, Tf 

and >* respectively. 


The Approach of Fenech and Rohsenow: 

Although the work of Fenech and Rohsenow [85] came later 
than the other two which remain to be discussed, it is pre- 
sented here because of the similarity of approach to the pa- 
per described above. 

Fenech and Rohsenow [85] , using the contact model de- 
picted in Figure 3, obtained solutions for the temperature 
distribution in the metals in region of the contact. How- 
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ever, since the exact boundary conditions are extremely com- 
plex, the solutions were found for a set of approximate or 
average boundary conditions. The resulting expression for 
the total contact conductance, i.e., the sum of the fluid 
and solid conduction terms is, 


$1 

S/I 




0-0 1 


(t I;)] 

1 -it, Ht, J 



(H-15) 


In which the quantities not previously defined are n = N/A = 
the number of contact spots per unit apparent area, ^ ~ 

the square -root of the ratio of the real to apparent contact 
areas, and the function f ( 6 ) is a transcendental func- 
tion. The value of f (€) is plotted in the paper, how- 
ever, the authors point out that for most practical situa- 
tions €<.0.1 and that in this range f^ 1.0. In equa- 
tion (11—15) it is noted that h is the sum of two frac- 
tions. The first fraction, shown by the square brackets in 
the numerator, represents the gap heat flow, and the second 


fraction represents the solid metallic conduction. In order 
to use equation (11-15) three remaining quantities must be 
evaluated, namely, $ M and € . The authors develop the 
following approximate expression, valid for small 6 , for 
the average fluid gap thicknesses. 
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S- = - ^ 7 " )C = /,Z ( 11 - 16 ) 

Where the two yU± are the volume average void thicknesses of 
the two surfaces. The parameters yu , 7t and € may be 
determined from surface profile measurements obtained with 
a prof ilometer . For each surface it is necessary to obtain 
two profiles taken perpendicular to one another. If the two 
surfaces are not randomly rough the orientation of the two 
profile readings on the two surfaces should correspond to 
the orientation they will have when contact is made. The 
recorded profiles are reproduced on transparent sheets, and 
all measurements are made visually and/or graphically by 
superpositioning the corresponding profiles. The volume av- 
erage void thickness is determined with the aid of a planim- 
eter. An actual count is made to determine n. Actual con- 
tact area, A , is measured directly, thus £ can be calcu- 
c 

lated. The profiles are shifted laterally a small amount 
and the processes repeated. Averaging the results thus ob- 
tained gives better values for the measured quantities. 

Test results and comparison with the above approach are dis- 
cussed in the next section. 
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The Approach of Laming: 

The same basic contact model as is shown in Figure 3 
was employed by Laming [l43] in his simplified approach to 
predicting solid conductance. Laming assumed, as did Cetin- 
kale (Veziroglu) and Fishenden, that when a contact is load- 
ed the softer metal yields plastically until the local pres- 
sure borne by an a-spot is equal to the Meyer hardness val- 

ue, M, of the material. 

Thus from the simple force balance, 

PA = MA C » M (Wn-cp) 

a relation for the spot radius, a, can be found. 

/i. _P\i 

4 ~ (firr m J (n-17) 


In equation (11-17) , n is the a-spot density, or number of 
contacts per unit apparent area (N/A) . Laming then used the 
work of Holm (l 18) who showed that the ideal constriction con- 
ductance of a single a-spot is 2ak , i.e., for a sincrle 

m 3 

a-spot which is infinitesimally small compared to the heat 
channel feeding it and is adjacent to a non-conducting gap. 

, the author assumes that this relation can be modified 
to apply to a single finite size a-spot by use of a "con- 
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The term (1-f) , which is dimensionless, is actually the 
first two terms of a series derived by Roess 0-84] . How- 
ever, the series is a power series in (a/b) having only odd 
terms, and, as Laming points out for most practical cases 
only the first two terms are needed. Thus, from Roess ‘s 
work f has the approximate value 1.41 (a/b), see equation 
(11-22) . In order to evaluate n. Laming assumes that the 
waviness of each surface consists of parallel ridges (such 
as would result from using a shaper) of wavelengths X and 
X 2 > an( ^ that the angle of intersection between the ridge 
lines of the two surfaces when they are brought together is 
o(. . Then the conductance is given by 

/ / -ain* P \ x 

\= hT (oTX'M j f 11 - 18 ) 

For the fluid gap conductance Laming used k / $ , where 

f f 

is the effective fluid gap thickness, and its value must 
be determined from experimental data. In his work Laming 
found the value of $ to be equal to two-thirds of the vol- 
ume average gap thickness, which can be found from surface 
profile records. His method involved determining the zero 
load conductance, by extrapolation, for tests of the same 
set of surfaces with fluids having different thermal conduc- 
tivities. Only the parameter (a/b) remains to be found. By 


forming the ratio of solid conduction heat flow to total 


heat flow. Laming shows that the ratio (a/b) is given by 

a = (P d* ) * 

«s I* T,J 

Thus the final form of the total contact conductance is 


^ -f- , P 1 1 

i/ /-/ ( tt A, A» M J 


(11-19) 


where h s is given by equation (11-18) , and the 
ing there is given by 


-f= 1.4-1 



appear - 


(11-20) 


From which it can be seen that an iterative process is re- 
quired to find h since the equation for f contains h . 

s 

Laming points out that f is small and convergence is rap- 
id for most practical cases. Therefore, using Laming 's sim- 
piifisd approach, one needs only to have the Meyer hardness 
(approximately 3 times the compressive yield strength if un- 
available) and surface profiles which give both roughness 
and waviness, in addition to the thermal properties, to be 
able to predict contact conductance. 


The Approach of Clausing and Chao: 

Clausing and Chao ^58]j made analyses of both the mi- 
croscopic and macroscopic constriction conductances. For tha 
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microscopic constriction conductance analysis they employed 
the model shown in Figure 3 and used the work of Holm [ll8] 
and Roess [l84j , as Laming [143] did, to arrive at the fol- 
lowing relation for h 

ui-21) 

S *5 ( /b ) 

In which all terms have been previously defined except 
g (a/b) , which is Roess's series: 

= I - [.HOIIS (*A) + 0.21*71 ( V6)*+ • * • (II-22) 

The authors assume, following Holm [118] , that the a-spots 
do not deform completely plastically. They argue that the 
asperity deformation is partially elastic. To account for 
this they assumed that the average pressure borne by an 
a-spot is only a fraction of the microhardness, M. That is, 

P = i W , (o£ t, ^ I) (11-23) 

Finally, the authors assume that an average value of % = 0.3 
and a value of unity for g(a/b) are representative for 
most contacts. With these assumptions equations (11-21) and 
(11-23) give 

^ 3 P 

~ aj^r (II " 24) 
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Holm [l 18] demonstrated experimental justification for the 
assumption of equation (11-23) , however, he recommended a 
value of K — 0.5. The lower value assumed by Clausing and 
Chao is hard to justify. 

For their analysis of the macroscopic constriction re- 
sistance Clausing and Chao employ basically the same model 
as shown in Figure 3, except that here the contacting region 

radius is called a and the heat channel radius is called 

Li 

(L for large). The contacting portion (radius a ) con- 

sists of a large number of contact spots (each with a radius 

a) and there are no contacts outside the radius a . Figure 

L 

6 shows the geometry employed in the case of a single macro- 
scopic contact. Again employing the work of Holm and Roess 
the authors arrive at the following relation for the solid 
conductance : 


1 "I>l 


(11-25) 


For this macroscopic problem the authors assume that the 
contact region area is controlled by elastic deformation of 
the contact members. They assume that the flatness devia- 
tions of the surfaces may be simulated by spherical caps of 
radii Y^ and Y as shown in Figure 6. The heights, d 
and d 2 , of the unloaded caps are called the "equivalent 



Microscopic Contacts 



Macroscopic Constriction Model 
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flatness deviations” of the respective surfaces. To obtain 
the load-area relationship for this model the authors use 
the well known elastic deformation work of Hertz £L16j . 

This results in the following relationship for ( aL /bL) : 

= !. 2 Ss [ (11-26) 


In which E is the harmonic mean of the two moduli of 
m 

elasticity, and E^. The dimensionless quantity in 

brackets in equation (11-26) is called the "elastic conform- 
ity modulus" and given the symbol, ^ . Thus, substituting 
(11-26) into (11-25) gives 


4 ^ = , 4 (x) 

rr C/.a»s*?v 


(H-27) 


It should be noted that b was a known quantity in the 

L 

work of Clausing and Chao, since their test specimens were 
made to approximate the model shown in Figure 6. For any 
practical contact surface b^ would be difficult to evalu- 
ate. The present writer assumes that b^ could possibly be 
estimated from surface waviness. The authors compare the 
magnitudes of the microscopic and macroscopic constriction 
resistances by calculating values of the ratio of equations 
(11-27) and (11-24) for some "representative" surfaces and 
pressures. From these calculations they conclude that for 



39 


many surfaces the macroscopic constriction resistance is 
dominant by approximately 2 orders of magnitude. The re- 
sults of a comparison of their theory and experiments is 
discussed in section III. 

Summary on Solid Conduction Mode 

An evaluation of the relative merits of using any of 
the above four methods hinges on three things (1) type of 
surface; (2) type of data available; (3) accuracy de- 
sired. 

For reasonably flat, rough surfaces the method of 
Cetinkale (Veziroglu) and Fishenden, Fenech and Rohsenow, 
and Laming are all applicable. Of these three Laming’ s ap- 
proach is the simplest and requires the least amount of data 
and effort. The method of Fenech and Rohsenow requires much 
more work in evaluating some of the parameters but gives 
more accurate results. The Cetinkale (Veziroglu) and Fish- 
enden approach probably lies somewhere between the other two 
in accuracy and has the disadvantage of requiring some ex- 
perimental heat transfer results, unless the average values 
of the parameters which they report are used. Clausing and 
Chao jj>8, pa. 59] state that ".....it is doubtful that their 
models [[models which consider only microscopic constriction] 
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can be valid." However, the above three approaches have met 
with reasonable success in predicting contact conductance of 
surfaces for which they are applicable. 

For smooth surfaces in which there is considerable 
large scale waviness there is no doubt that the emphasis of 
Clausing and Chao on the macroscopic constriction is appli- 
cable. Their experimental results for cylindrical specimens 
with spherical ends agrees well with their theory. However, 
there is no indication of how their method could be applied 
to practical surfaces. For surfaces which fall in between 
the "rough-flat" and "smooth-wavy" categories no theory has 
yet been published. 


III. LITERATURE SURVEY 


The purpose of this section is to provide a survey of 
the current status of the contact conductance literature. 

It was the decision of the writer not to present the stand- 
ard type of review of this literature. Instead, it is 
hoped, that a more useful document would result from a pres 
entation based on the following outline: 1) reviews of 

some of the major works in the field and some of the high- 
lights of past research, and; 2) a tabulated review, cate 
gorized by the type of information to be found in the var- 
ious references. 

Discussion 

Most of the earliest published works concerning con- 
tact resistance were written by people who were interested 
in measuring thermal conductivity. In the course of their 
experiments these investigators were aware of the tempera- 
ture discontinuities caused by contact resistance at the in 
ter faces. For the most part these early writers were con- 
cerned with attempts to eliminate the contact resistance. 

The earliest work found by the present writer that 
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deals with thermal contact resistance was published by 
Northrup 0-71] in 1913. Northrup’s paper is an exposition 
of the analogy between the flow of heat and the flow of 
electricity, which leads up to a proposal of an apparatus 
for measuring thermal conductivity or thermal contact re- 
sistance. The method is based on comparing the temperature 
drop across a fixed length of a standard material to the 
temperature drop across the same length which consists part- 
ly of the standard material and partly of an unknown mate- 
rial or a contact resistance. The results are expressed as 
an equivalent length of the standard material. Northrup re- 
ports the contact resistance of a single copper-copper in- 
terface as being equivalent to 31.2 cm of copper at a load 
of 1.6 kg on the 3.8 cm diameter rods. Although Northrup’s 
method left a lot to be desired one could not find ther- 

mal conductivity accurately unless the contact resistance 
was eliminated his idea of expressing thermal contact re- 

sistance as an equivalent length of some material was to be- 
come standard practice for a while. 

In 1919 Taylor C2173 was interested in measuring ther- 
mal conductivities of some insulating materials for electric 
motors. Taylor, by experiment, realized the shortcomings of 
Northrup's relative method and modified the apparatus to 
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eliminate the problems. The method used by Taylor is a com- 
parative method in which the temperature drop in each part 
of the apparatus is measured separately. It is basically 
the method which has been used by almost all thermal con- 
tact resistance investigators since then, including the 
present work. Taylor, does not report any values of contact 
resistance but, does report that the addition of such sub- 
stances as glycerin, vaseline, glue and shellac reduced the 
contact temperature drop to a negligible amount in the cases 
of low heat flux. 

At the Bureau of Standards, in 1922, Van Dusen [ 225 ] 
developed the basic apparatus one further step to include a 
spring for controlling the pressure of contact between two 
or more thermal conductivity specimens in series. Van Dusen 
was concerned with measuring thermal conductivities of met- 
als and very thin slices of poor conductors. Thus his in- 
terest in contact resistance was in finding ways to elim- 
inate it, and he reports that wetting the surfaces with wa- 
ter or light mineral oil effectively eliminated the prob- 
lems. He states that for his standard specimens (brass) 
with surfaces ground fairly flat, actually, convex with a 
of curvature of 30 meters, the dry contact resist- 
ance was equivalent to 1 - 1.5 cm of brass at contact pres- 
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sures of 1-5 kg/cm 2 . He also reports that whan the sur- 
faces were lapped together with water or dilute glycerin the 
contact resistance was reduced to the equivalent of less 
than 0.5 mm of brass. Although the data are scant it is in- 
teresting to note that Van Dusen was aware that surface fin- 
ish was also important to contact conductance. 

Jacobs and Starr [127] , in 1939, published the first 

work in which the sole purpose was to investigate thermal 
contact resistance. It is also, oddly enough, the first 
work for contacts in a vacuum environment. They were inter- 
ested in obtaining data for designing a mechanically oper- 
ated thermal switch for use in low temperature research. 
Their data, which were obtained for gold, silver, and copper 
surfaces polished to approximately optical flatness, were 
reported as curves of contact conductance versus contact 
pressure. In their apparatus the contact pressure was var- 
ied by controlling the pressure in a bellows which was ex- 
ternal to the vacuum system. Measurements were made at room 
temperature (25°C) and at - 195°C and for contact pressures 
from zero to 2.5 kg/cm . The most interesting result of 
this work was that they observed a linear variation of con- 
tact conductance versus contact pressure for copper, but for 
silver and gold the variation was approximately with contact 
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pressure to the 1/3 power. Jacobs and Starr were obviously 
aware of the primary variables affecting contact resistance, 
however, they failed to report any actual, measured surface 
conditions or vacuum pressure. 

In 1948, Fishenden and Kepinski [86] published a short 
note in which they reported the contact resistance of an in- 
terface formed by bringing together the two surfaces of a 
round steel rod which had been cut through with a saw. The 
contact resistance is given as being equivalent to approxi- 
mately a 0.001 inch thick layer of air when the rods were 
replaced together in the orientation in which they were 
sawed, and a .002 inch thick air gap when the two pieces 
were rotated 30° about the axis of the rods. No information 
is given about the steel composition, surface roughness or 
contact pressure. Although the work was called a "prelim- 
inary note" by its authors the present writer has been un- 
able to find any subsequent published work to follow up the 
original note. 

The years 1948-49 represent the beginning of the sig- 
nificant work in contact resistance. Interest in the subject 
began to increase, and the caliber of the work began to im- 
prove at that time. 

Brunot and Buckland [47] experimentally investigated 
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the contact resistance for two types of steel contacts. Re- 
sults are given for 2 -inch square laminated steel blocks 
(with heat flow parallel to the lamination planes) for tests 
in which the lamination planes of the two blocks were par- 
1 snd perpendicular but in direct contact, and for tests 
in which thin shims of aluminum, steel or cement was placed 
in the contact. Results are also presented for two 2-inch 
square cold-rolled steel blocks with surface roughnesses of 
4 to 1000 microinches (rms) . Contact pressures in the first 
case ranged from 25 to 200 psi and in the second case the 
range was from near zero to 300 psi. For the laminated 
blocks the authors found that the parallel laminations re- 
sulted in higher resistances than the perpendicular lamina- 
tions. It was also observed that the addition of the thin 
shims served to reduce the contact resistance — with alumi- 
num foil giving the greatest reduction. For the solid steel 
blocks it was found that the smoother the surface finish the 
lower the contact resistance, and that the variation of con- 
tact resistance with pressure is greater for the rougher sui^ 
face finishes. 

Weills and Ryder [233] presented the most complete set 
of data that had been published at that time. They present- 
ed contact conductance results for dry (air in the interface) 
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and oil-filled interfaces between surfaces of various rough- 
nesses of aluminum, steel, and bronze. Test blocks were 3 
inch diameter bars 3 inches long. Heat flow was measured 
independently and thus allowed the thermal conductivities of 
the specimens to be measured. The contact conductance re- 
sults show the effects of mean joint temperature, contact 
pressure and surface finish. The effect of plating the 
steel surface with copper was also studied. Complete infor- 
mation about the test samples, including chemical composi- 
tion and mechanical properties, and test conditions are giv- 
en. From their results the authors found that the thermal 
conductance of a dry joint increases with contact pressure, 
linearly for steel, and exponentially for aluminum and 
bronze, and that the thermal contact conductance increased 
with decreasing surface roughness for both dry and oil-filled 
contacts. At low contact pressures the thermal contact re- 
sistance was decreased by a factor of 2 with the addition of 
oil in the interface. They also reported that the contact 
conductance increases slightly with mean interface tempera- 
ture, the increase being greater for smoother surfaces. 
Weills and Ryder are also apparently the first to report a 
hysteresis-like variation in the contact conductance versus 
contact pressure relation when the contact pressure is de- 
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creased following an increase. The decreasing pressure con- 
ductances were found to be higher than those measured at the 
same pressure when reached by an increase in pressure. 

The effects of contact pressure, mean interface tem- 
perature and surface roughness of a series of steel— steel 
joints were investigated by Kouwenhoven and Potter [1403 . 

In this paper thermal contact resistance was measured for a 
contact formed between two specimens of mild steel, one 
specimen in every test having a surface roughness of 3 mi- 
croinches ( rms) while the second specimen surface roughness 
was varied over the range of 3 to 3320 microinches (rms). 
Tests were made for several combinations of roughness with 
contact pressure variations from 195 to 3000 psi while keep- 
ing the mean joint temperature fairly constant (170 - 197°F). 
A similar set of tests were made for a contact mean joint 
temperature of 600°F. The former set of tests was made with 
air in the joint, the latter set was done with Argon to pre- 
vent corrosion. Two other sets of contact resistance data 
were obtained for several combinations of surface roughness 
by varying the heat flux so that the mean joint temperature 
was varied from 350 to 700°F while keeping the contact pres- 
sure constant at 195 psi in one case and 1575 psi in the 
other. The authors concluded that the thermal contact re- 
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sistance decreases exponentially with contact pressure .with 
the rate of decrease being greater for rougher surfaces. 

For very smooth surfaces they found that the contact resist- 
ance was essentially independent of pressure. Both of these 
conclusions are in qualitative agreement with other investi- 
gators [24,47,233]. They concluded that the contact resist- 
ance at a constant contact pressure is substantially inde- 
pendent of contact temperature (over the range investigated). 
However, their data contradict this in several cases in 
which the variation was about the same as that found by 
others. The authors also have an interesting comparison be- 
tween a theoretical contact area ratio and the actual ther- 
mal resistance ratio for changes in contact pressure in the 
case of some specimens with ruled roughness. For the simple 
theoretical model the ruled ridges are assumed to be 45° 
isosceles trapezoids in profile. The comparison for a 3 mi- 
croinch surface was very good, but for the rougher surfaces 
the thermal resistance ratio decreased more rapidly (with 
increasing contact pressure) than the contact area ratio. 

One of the more significant theoretical works in the 
field was published in 1951 by Cetinkale (Veziroglu) and 
Fishenden [55] . The details of their approach were discussed 
previously in section II. The authors presented a drawing 
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of the test apparatus but do not present their test data. 
They mention only the ranges of the test parameters and in- 
dicate that agreement between theory and experiment was good. 

Barzelay, et al [24 - 29] published a series of NACA 
reports on their investigations of thermal contact conduct- 
ance. Their results, too numerous to detail here, covered 
both cut-bar apparatus measurements and conventional skin- 
stringer type aircraft joints. Materials tested were alumi- 
num (7075-T6 and 2024-T4) and 416 stainless steel. A large 
range of surface roughness and contact pressure was studied. 
The effects of mean joint temperature and the addition of 
thin foils to the interface were studied. Complete test 
conditions and results are reported. The most interesting 
results of these references are the following. They are ap- 
parently the first to observe that surface roughness is not 
always dominant in controlling contact conductance but that 
flatness deviation (waviness) sometimes predominates. It 
was also noted that in some of the skin-stringer tests the 
warping of the members due to thermal expansion produced 
some adverse effects on the contact heat transfer. 

An equally important "first" reported by Barzelay, et 
al [25] was the difference observed in contact conductance 
for dissimilar metal joints when the direction of heat flow 
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was reversed. They found that the contact conductance for 
an interface in which heat flowed from an aluminum specimen 
to a stainless steel specimen was several times higher than 
when the heat flow direction was reversed, other test con- 
ditions being the same. 

Barzelay and Holloway ([26,27] are the first published 
work to study the effects of an interface on thermal trans- 
ients. They made tests on a number of riveted aluminum 
stringer -skin combinations in which they measured the tem- 
perature at several points on a T-shaped stringer riveted 
to a flat skin. The skin was subjected to a constant ra- 
diant heat flux to simulate aerodynamic heating. It was 
found that the interface resistance had a significant de- 
laying effect on the time-temperature history of a given 
point on the stringer. There was considerable scatter in 
the conductance data which they attributed to the warping 
effects. 

Further work on the directional effect in dissimilar 
metals was carried out by Rogers [187] . In Rogers' experi- 
mental apparatus the direction of heat flow was changed with- 
out separating the test specimens. Tests were carried out 
with aluminum and steel specimens at a constant contact pres- 
sure of 122 psi. Results showed that for the same aluminum 
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steel joint (in air) the contact conductance was about 20% 
higher when heat flowed from the aluminum to the steel than 
when the heat flow was reversed. When these specimens were 
tested in a vacuum the difference was approximately 100% 
due to the low conductance in the vacuum. Other combina- 
tions tested were chrome 1 -a lumel and copper-steel. The for- 
mer showed no clearly defined directional effect, the latter 
only a slight effect. Rogers suggested that the directional 
effect might be associated with the mechanism of conduction 
at the points of metallic contact. 

Williams [236] , in a comment on Rogers' paper, sug- 
gested that the directional effect was a result of surface 
oxidation. Williams states that the difference in lateral 
(parallel to contact plane) thermal expansion could liter- 
ally "scour" the oxide layer off the aluminum surface. In 
response to Rogers’ paper. Moon and Keeler [167] applied the 
theory of solid state heat conduction. They showed that the 
effect could be qualitatively explained by demonstrating 
that the electronic conduction contribution was directional. 
They stated that accurate quantitative-analysis was impos- 
sible due to lack of data on the work functions of metals. 
However, Ashby Q>] claims that Moon and Keeler made an error in 
their analysis' and that such an effect is not demonstrated 
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by the corrected calculations. It is clear that at this 
time the directional effect in dissimilar metal contacts is 
not well understood. Results of research which is presently 
in progress at several institutions may serve to clarify 
this phenomenon. 

In the design of nuclear reactors contact resistance 
can be an important consideration. Several papers have been 
written which are concerned with the contact resistance of 
materials and/or configurations which are applicable to re- 
actor design. One of the first of these was published by 
Boeschoten and Van Der Held [l4lj in which surfaces of alu- 
minum-aluminum, aluminum-steel and aluminum-uranium were 
studied. One of their purposes was to determine the gas and 
solid conduction contributions separately. This was accom- 
plished by running tests at constant contact pressure and 
mean interface temperature while varying the ambient gas and 
gas pressures. The conductance curves for three gases, air, 
helium and hydrogen, were extrapolated to zero gas pressure, 
thus giving the value of the solid-solid conductance. From 
a simplified analysis they calculated the average a-spot ra- 
dius to be approximately 30 microns. This value was found to 
be relatively independent of materials and contact pressure 
for the range of their test conditions. 
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Wheeler |J2 353 wr °te a survey paper in which he exam- 
ined some of the published vacuum conductance data. He con- 
cluded that most of the existing (unclassified) work could 
not be extrapolated with sufficient accuracy to the higher 
heat flux levels characteristic of nuclear reactors. In a 
later work £234] Wheeler presented some data for materials 
and heat flux levels typical of nuclear reactors. He also 
tried a different approach to conductance data correlation 
which consisted of plotting contact conductance versus the 
ratio of apparent contact pressure to the yield strength of 
the softer contact material. His results showed consider- 
able scatter. He also attempted to find the effective gap 
thickness by conducting tests on the same joint in a vacuum 
and with different gases. The results were highly doubtful 
since the gap thicknesses calculated were on the order of 2 
to 10 times the total surface roughness heights. It seems 
likely that these results may indicate the presence of large 
scale flatness deviations of which Wheeler was unaware. 

Skipper and Wootton £203] studied the contact conduct- 
ance of uranium-Magnox and uranium-aluminum joints. In ad- 
dition to the effects of mean temperature and gas pressure 
they also reported data on the effects of very thick oxide 
films. For gas-cooled reactor applications Sanderson £l8EQ 


also investigated the contact conductance between uranium 
and Magnox surfaces. Data are reported for the effects of 
contact pressure, interface temperature, ambient gas pres- 
sure and surface finish. The reduction of contact conduct- 
ance by oxide and nitride surface layers is also reported. 

One of the first attempts to make a general correla- 
tion of thermal contact conductance data was made by Graff 
p.03j . He plotted the data of several investigators in the 
form of two dimensionless groups: hp/k ^ versus p/B. The 

former group is the product of conductance and contact pres- 
sure divided by the product of the thermal conductivity and 
the density of the metal. The latter group is the ratio of 
contact pressure to Brinell hardness number. These dimen- 
sionless groups did not correlate the data. Graff suggested 
the reason for this was the lack of a roughness parameter. 

Most contact conductance data published were obtained 
with cylindrical specimens with axial heat flow. The work 
of Barzelay [24-29] on "practical" joints has already been 
mentioned. Aron and Colombo [13] report some data on a sin- 
gle, bolted type aircraft joint. One of the first experi- 
mental works using flat plate specimens was published by 
Fried and Costello [96] . They measured thermal contact con- 
ductance for 5 x 5 x 1/8 inch thick plates of aluminum and 
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magnesium. Contact pressures ranged from 2 to 35 psi; a 
range the authors felt representative of bare space vehicle 
joints. They found that smoother surface finishes gave 
higher conductances, but that flatness deviations could have 
significant effects on the conductance for smooth and rough 
surface finishes. The authors also reported that the addi- 
tion of soft shim materials, such as aluminum and lead foils 
and thin, copper wire -mesh cloth, could improve the conduct- 
ance by as much as a factor of 2 - 4 times the value for a 
bare joint. Of these three, lead foil gave the largest in- 
crease. Fried, in a later paper £94]] , found that the ad- 
dition of a silicone grease could greatly increase the con- 
tact conductance . 

Jansson 0-28) also investigated interstitial filler 
materials for contact conductance improvement . He found 
that for the materials and conditions tested, indium foil 
gave the largest improvement. Other materials tested were, 
in order of descending improvement, epoxy cement, lead, alu- 
minum and gold foils. Indium foil was also found to give 
better conductance improvement than lead, and aluminum foils 
by Koh and John [j.393 . From the results of this study they 
concluded that foil softness was more significant than the 
thermal conductivity of the foil. They also demonstrated 
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experimentally that there is an optimum foil thickness for 
a given joint. 

Further work on thin plates was done by Stubstad 
C206,207j . He measured the contact conductance for a joint 
between two 3 x 3 x 1/8 inch plates in a vacuum at very low 
contact pressures, i.e., 2 to 20 psi. Materials tested were 
aluminum, copper and stainless steel. He observed that the 
contact conductances was an order of magnitude higher at at- 
mospheric pressure than in a vacuum of 10“^ mm Hg, which 
demonstrates the large gas conduction contribution at low 
contact pressures. It was also observed that test repeti- 
tion with different sample orientations produced large vari- 
ations in the contact conductance. 

In 1961, Laming [143] presented a simplified analyti- 
cal approach to the prediction of contact conductance. The 
details of his approach were presented in section II. He 
reported test data for steel-brass, steel-aluminum, brass- 
aluminum and brass-brass contacts with air, water and glyc- 
erol in the joints. Agreement of his test data with his 
theoretical equations was poor at low contact pressures 
(20 psi) and improved somewhat with higher contact pressures. 
However, he found that if a load dependent value of the 
Meyer hardness was assumed, all his test results correlated 
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with the theory very well. The load dependent microhardness 
hypothesis results in high hardness values at low loads. Al- 
though there is some evidence to support this hypothesis. 
Laming [143, p. 75] states that the only claim made for it is 

" its value in correlating the heat transfer data." 

However, his assumptions have at least been partially justi- 
fied by the recent results of Williams [238] . Williams ran 
tests on mild steel and nickel specimens and controlled the 
number of contacts spots by using ridged surfaces of various 
ridge frequencies. His test results showed that the effec- 
tive hardness of the nickel specimens increased by a factor 
of 5 when the average spot load was reduced from 1000 lbf to 
0.1 lbf. Similar results for mild steel were also reported. 

Laming' s single theory (i.e., without load dependent 
hardness) predicts that the solid conductance should vary as 
contact pressure to the one -half power. His test results 
show that the variation is nearer the two-thirds power. 
However, Fried £94] found, using his own data for aluminum 
and magnesium, that the contact conductance in a vacuum 
(therefore only solid conductance) varied very nearly as the 
one-half power of contact pressure. This may suggest that 
Laming' s discrepancies arose from his method of estimating 


the fluid conduction contribution. 
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Fenech and Rohsenhow [853 presented a thorough analyt- 
ical treatment of contact conductance, the details of which 
have already been discussed. To test their theory they ob- 
tained contact conductance data for an armco iron-aluminum 
joint with maximum roughnesses of 150 microinches (rms) . 
Contact pressures ranged from 92 to 2625 psi. Agreement be- 
tween measured results and theory was very good. They also 
tested several "idealized" contacts. These consisted of 
stainless steel specimens with a machined pyramid surface 
against an optically flat surface, and an iron pyramid sur- 
face against and an optically flat-topped pyramid surface 
of aluminum. Plots of measured and theoretical contact con- 
ductance versus mean joint temperature for contact pressures 
of 92 to 6226 psi $ho>w good agreement for these surfaces, 
the agreement being better for the stainless steel speci- 
mens. Another set, of tests are reported for some solid cyl- 
inders with a neck machined in them to simulate a single 
conduction channel and a-spot (see Figure 3). Results are 
plotted for three sizes of neck, i.e., three values of the 
ratio of spot radios to heat channel radius ( a /b) . Curves 
of conductance versus mean contact temperature for air, wa- 
ter and mercury in the gap show reasonable agreement between 


theory and experiment. 
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Henry [115] compared the test data of Adamantiades [ 4 ] 
with the theoretical prediction method of Fenech and 
Rohsenow. The experimental results were for stainless steel 
surfaces which had been ground to a "mirror finish" and sub- 
sequently blasted with small glass spheres to achieve ran- 
dom roughness. Plots of measured and theoretical contact 
conductance versus contact pressure for three mean tempera- 
tures were presented. Agreement was very good. 

The above comparisons show that for test conditions 
which meet the assumptions of the theoretical analysis the 
agreement with experiment is good. However, to date, the 
theory has not been compared with "everyday" engineering 
surfaces, mostly due to lack of sufficient test information. 
Thus there is still some question of its general applicabil- 
ity . The only real disadvantage to the approach of Fenech 
and Rohsenow [85] is the very involved graphical procedure 
that is required to determine the necessary surface parame- 
ters. Others at Massachusetts Institute of Technology have 
been working on methods of simplifying the work required t.o 
obtain these parameters. The results of these efforts are 
summarized by Henry [115] . 

Clausing and Chao [58} presented a theoretical analy- 
sis which assumes that macroscopic heat flow constrictions 
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caused by the flatness deviations of a surface dominate the 
metal contact resistance. This is opposed to the analyses 
of others [55,85,143] which assumed only microscopic con- 
strictions due to uniformly distributed microscopic contact 
spots. Clausing and Chao conducted experiments on cylindri- 
cal specimens whose contact surfaces were spherical caps. 
These specimens matched their theoretical model, which is 
shown in Figure 6. Materials tested were aluminum, stain- 
less steel, and brass. Contact pressures varied from about 
10 to almost 1000 psi. Agreement between theory and experi- 
ment was very good for all reported data. It should be 
noted that the test specimen and theoretical model have only 
one macroscopic contact. Most real surfaces of practical 
importance would have several such contact areas. The au- 
thors do not suggest how the theory could be applied to a 
practical surface. The real value of their approach, and 
their experimental verification, is that it does account for 
macroscopic effects, which are known to be important under 
some conditions. However, as is true for the other theoret- 
ical approaches, further work remains to be done before the 
contact conductance of general engineering surface can be 
predicted. 

In addition to the work of Barzelay [253 * which has 
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already been mentioned, the only other published experimen- 
tal work on thermal transient effects on contact conduct- 
ance is the recent paper by Schauer and Giedt [190] . The 
authors derived theoretical equations and devised an experi- 
mental method for determining the contact conductance be- 
tween two thin plates during transient heat transfer. The 
method is based on the heating of one of the plates with a 
capacitor -bank discharge and recording the temperatures of 
the contact surfaces as they come to thermal equilibrium. 
Tests were performed on aluminum-stainless steel and stain- 
less steel-ceramic interfaces. The heating time was approx- 
imately 100 microseconds and the temperature data were re- 
corded by an oscilloscope for about 160 milliseconds. The 
metal specimens were 0.032 inch thick and the ceramic speci- 
men was 0.302 inch thick. These test results showed that 
the contact conductance of the aluminum-stainless specimens 
increased sharply with time, whereas the opposite was true 
for the stainless steel-ceramic specimens. The method of 
calculation assumes that no heat is transferred during the. 
short capacitor discharge time, that all the heating occurs 
in the plate of lowest electrical resistivity and that an 
instantaneous temperature rise was produced. The first two 
of these assumptions seem to have been satisfied by the ex- 
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periment. However, due to the fact that the thermocouples 
were embedded in epoxy cement there was an indicated tran- 
sient rise in the temperature difference. This made it nec- 
essary to extrapolate to time zero to get the initial tem- 
perature difference, and also ruled out experimental verifi- 
cation of a doubtful boundary condition. This together with 
possible errors in thermocouple location make the results 
highly doubtful. Even if the results are valid they indi- 
cate that the conductance approached a steady value very 
rapidly (say 100 milliseconds) and therefore the indicated 
transients would be of no significance in most practical 
situations. The method also has the strict limitation of 
being applicable only to very thin pairs of dissimilar sol- 
ids due to the method of heating. 

One other work on transient conductance effects de- 
serves mention here. Aaron and Blum { 2 } performed a theo- 
retical analysis of the effects of varying the ambient gas 
pressure on the contact conductance and temperature distri- 
bution in two contacting cylinders. They predicted the ex- 
istence of a threshold pressure above which the contact con- 
ductance would be independent of pressure. They also demon- 
strated that for most practical joints the gas pressure with- 
in the contact voids would respond very rapidly to ambient 


64 


pressure changes. Some results of the present work which 
appear to be at least a partial experimental verification 
of these predictions will be discussed in section VI . 

Summary 

A glance at the above brief review shows that the phe- 
nomenon of thermal contact resistance is far from being com- 
pletely understood. Although a great deal of experimental 
data have been obtained, empirical correlations of them have 
been unsuccessful. Theories have been presented, and test 
data on well prepared samples of particular configurations 
have agreed well with the theories. Yet, no means is pres- 
ently available for the accurate prediction of the thermal 
contact conductance of a general engineering interface. The 
failure of past work is due mainly to the fact that no sta- 
tistically meaningful means of characterizing the contact 
between two surfaces is yet available. The validity of the 
statement that the past work is truly a failure is supported 
by the fact that the National Aeronautics and Space Adminis- 
tration has recently awarded a contract to E. Fried of Gen- 
eral Electric to measure, individually, the contact conduct- 
ance of over 100 separate practical joints for the Apollo 


spacecraft. 


65 


Tabulated Review 

The brief reviews presented above obviously do not 
represent more than a fraction of the total number of ther- 
mal contact resistance references. However, they do repre- 
sent, in the writers opinion, a good sampling of the liter- 
ature, both from the standpoint of historical development 
and the standpoint of indicating the important effects that 
have been observed. An interested reader can refer to the 
tabulated review below to obtain more information about any 
of the particular aspects of contact resistance. 

There are 244 references listed in the Bibliography 
section which pertain to thermal contact resistance or a 
closely related topic. No claim of completeness is made. 
However, it is believed that at least a majority of the sub- 
ject references are included in this list. The results of 
the review are presented in Tables 1,2 and 3. The headings 
used in the three tables are discussed individually below. 
Two points concerning the tables, which the writer wishes 
to emphasize are: 

1. The categories were chosen which seemed to be the 
most useful based on the writer's own experience. 

2. The writer was unable to obtain copies of some of 
the references listed in the Bibliography. There- 
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fore some references do not appear in any of the 
tables. Only those references which were reviewed 
by the writer are included in the tables, with a 
few exceptions where the title clearly indicated 
one of the special or related topics. 

Discussion of Reference Tables 

Those references which contain experimental results 
for contacts in the presence of a conducting fluid are 
listed in Table 1. References which contain data on the 
specific subheadings are listed next to each subheading. 

The numbers in all the tables refer to the corresponding en- 
tries in the Bibliography. The subheading topics are self 
explanatory. 

Table 2 contains those references which present exper- 
imental results for contacts at low ambient pressures. Most 
of these are, of course, for air vacuum. Subheading topics 
are the same as for Table 1. If a reader is interested in a 
combination of topics, e.g. bolted and riveted joints with 
interface fillers, it is a simple process to check both sub- 
headings and find the common entries. 

Table 3 contains references pertaining to special or 
related topics. The first subheading contains those refer- 
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ences which have anything that might be considered a the- 
oretical approach to contact resistance. References that 
allude in any way to thermal transients in the presence of 
thermal contacts are in the second subheading. The third 
subheading is a list of references which contain theoretical 
analyses of the temperature distribution in composite solids 
of various shapes and with various boundary conditions. Al- 
most all of these assume perfect thermal contact between re- 
gions, but are included here as a closely related special 
topic. The remaining related topic subheadings are self ex- 
planatory. 

As mentioned above, no claim of completeness is made 
for the entries in the tables. However, checking those ref- 
erences listed under any specific heading should provide a 
very good start. 


IV. THEORETICAL STUDIES 


A general description of the class of boundary value 
problems for which solutions are obtained in this work is 
presented first. This is followed by an outline of the der- 
ivations and the solutions obtained. The method of obtain- 
ing the numerical solutions for the same problems is dis- 
cussed briefly. Some remarks concerning the accuracy of the 
results obtained, including a comparison of the analytical 
and numerical results are given. A presentation and dis- 
cussion of the results obtained from a parametric computer 
study using the solutions concludes this section. 

The Boundary Value Problem 

The problem of interest here is that of obtaining the 
space -time temperature distribution in a composite solid 
medium consisting of two regions separated from each other 
by an interface which has a resistance to heat flow. 

Figure 7 depicts the geometry employed and some of the 
necessary nomenclature. Both regions are assumed one-dimen- 
sional and homogeneous within themselves, each with a con- 
stant thermal conductivity and thermal diffusivity. However, 


71 



72 



Fig . 7 . — Boundary Value Problem Geometry 
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there is a discontinuity in these thermal properties across 
the interface which necessitates separate solutions for the 
two distributions, as indicated in Figure 6 [ 248 ] . Both so- 
lutions must obey the one -dimensional Fourier equation. 

Thus for the distribution we have. 
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Applicable boundary conditions are as follows. Ini- 
tially both regions are assumed to have time -independent 
temperature distributions. To complete the problem descrip- 
tion it is necessary to specify the changes that occur at 
time 0 = 0 at any of the three physical boundaries and 

the other conditions imposed. For the present work the con- 
ditions at the outer boundaries (x = 0 and x = L) may be of 
the temperature or gradient type and the interface boundary 
will be of the gradient type. The detailed conditions em- 
ployed in the separate cases studied are developed below. 

Case A. 

Initially both regions are at a uniform temperature, 
t^. At & = 0 the temperature at x = 0 is suddenly raised 
to a value t_ and held constant. The temperature at x = L 
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is held at its initial value and the contact resistance at 
the interface is held constant. Symbolically, these condi- 
tions are as follows: 


±,(x,o) = i V 


(IV-2) 
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at the interface 
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The solution obtained from equation (IV-1) by the 
standard technique of separation of variables in product 
form is well known. With the proper choice of sign on the 
separation constants the solutions for the two regions are: 
Region (I) 

y yr * 8 /-*/* " ' (iv-9) 

~ c«, 

Region (II) 

= A x (L-x) + J F >^ S ^ X + 8 **J c **** ° 


(IV-10) 
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These solutions have been written as the sum of 
steady state and transient parts to take advantage of the 
fact that the steady state solutions are simply linear func- 
tions of x. For the transient portions the summations have 
already been indicated in anticipation of the requirements 
of expanding the initial boundary condition in a Fourier 
series. The constants A, and A„ are the well known 
steady state temperature slopes, and can be written as. 
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The solution of the problem is completed by determin- 
ing the unknown coefficients B_, C , F , and G , and 

n' n n n' 

the eigenvalues and $ n * It is noted that the form of 

coefficients was chosen in anticipation of simplifying the 
calculations. 

Application of the boundary condition (IV -4) leads to 
the following 

0= %.8-hCC* <') +°J 

Therefore, C = 0. 

n 

Upon applying (IV-5) it is seen that 
y% 
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G-h * -p^ W S^L 


(IV-13) 


Boundary condition (IV -8) is simply a statement of the 
assumption that there is no heat storage in the infinites- 
imal thickness of the interface. Its application yields 


S„ = Y„ (*•/«>)* <IV-14) 


Next, (IV-6) is applied and results in 
^ ^ ^ = ^ ^ ^t| ft^ 

' 7J 

Substituting (IV-13) and (IV-14) into the above equa- 
tion and applying some trigonometric identities gives 


F n = 


C++ £*£ 


(IV-15) 


The application of (IV-7) yields 

■A, |? s * y . K* = 

+ 3,^Ko-J 

Substituting (IV-13) , (IV-14) and (IV-15) into the 
above, together with some algebraic and trigonometric sim- 
plification, gives the following 

— ^ ^ (¥**) (IV-16) 


Equation (IV-16) is the eigenvalue equation from 
which the region I eigenvalues, n , may be found. Once 
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found, these values give the region II eigenvalues, § , 

n 

by means of equation (IV -14) . The coefficients F and 

n 

G n are then determined by equations (IV-15) and (IV-13) . 

Only the primary series coefficients B n remain to 
be found to complete the solution. It is well known that 
their values are determined by satisfying the so-called ini- 
tial condition, i.e., the time boundary condition. However, 
ordinary Fourier series expansion is inadequate here because 
the solution eigensets, namely, 

(IV-17) 

L 

are not orthogonal over the full x-interval (0, L) . How- 
ever * Tittle has shown [222] that these functions are "qua- 
si “orthogonal" and that an orthogonal set g n (x) can be 

constructed from a linear combination of the f (x) . These 

n 

orthogonal functions, g^ (x) , can then be used for the ex- 
pansion of the boundary function to determine the coeffi- 
cients B r . Proof of this theorem, called the "theorem of 
qua si -orthogonality, " is given by Tittle [222]. Briefly, 
the procedure is as follows. 

The set g^ (x) is given by 
^ -n tx ) - 


LX) * a.£x*L 


( IV-18) 
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In the above, the term C, is the "orthogonality factor." 
Application of the theorem results in 

D JT/' I+A,x)4m\ X C Z f^A 2 (*-*■)[%, Co4gj(JJx 

B n = 77, t (IV-19) 

X J* + c +G*c+4S n xfi/x 

The orthogonality factor, C, is determined from the 
orthogonality condition, which is, 

f t***)}**** 2*<*>J* = & *»*? r (IV -20) 

The weighting function, w (x) , appropriate to the coordi- 
nate system for this case is unity. Substitution of equa- 
tion (IV-18) into equation (IV-20) gives 

*x**Yt,xJx+- C 2 Jfa 4U, +G^C4*M]- ( iv-21) 

* *n* g h * J thf 

From this it can be seen that, alternatively, C 2 could be 

thought of as a discontinuous weighting function for the 

original functions f (x) . In a later and more general un- 

n 

published paper C 2 533 Tittle has given this interpretation. 
Regardless of interpretation equation (IV-21) will determine 
the value of C. However, Tittle has shown that the value 
of C is independent of boundary conditions and is the same 
for this class of problems, namely 

(IV-22) 
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The writer has verified equation (IV-22) by carrying 
out the integration of (IV -21) . With the value of C given 
by equation (IV-22) the integration of equation (IV-19) can 
be evaluated. 

The resulting form for the B is 

_ 2 n 
a* 

where D n is given by 




-jr 

Z "x 



2 




(IV-23) 


Assembling the results gives the final form of solu- 
tion for the distributions, which are. 

Region I 


T, \U.e)= \-A,x-£ 




(IV-24) 


Region II 

^ _ r ^ g 

T, (X i 9) a A x (L-X)~£ 2 <**$■**] € * ( IV-25) 

Wherein all quantities have been previously defined 
except the following 


T. 

Lp * t C 


ti <*■*> z 

-t. “ ’ti 


(IV-26) 


Equations (IV-24) and (IV-25) have been used to gen- 
erate theoretical data for some representative systems. 


80 


This data will be presented later in this section. 

Case B. 

Initially a steady-state temperature distribution is 
established through the system. At time 6 - 0 the contact 

conductance is suddenly changed to a new constant value, h^, 
while the outer boundary temperatures are kept constant at 


and t^. Symbolically these 

boundary conditions 

are. 


c cc 

(IV-27) 

TaC*,*)^ K (*•-*) 

X* L 

(IV-28) 

T t ( 0,9 \ = 1 

e>o 

(IV-29) 

T a (L,©) = O 

eyo 

(IV-30) 


x-ck, &>o 

( IV-31) 


&>0 

(IV-32) 

Boundary condition (IV-8) 

is also applicable. 

The 


method of solution is identical to that of the previous case 
and is not repeated here. Except for the numerator of the 
B n coefficients, the solution form is identical. These so- 
lutions are. 
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Region I 


T.cx.e)- I-*!* + 


Region II 


£ ^ ' ( V ~j0_£*± v 33) 


lia/ 


i ■» * * D_ 


(IV-34) 


Wherein F n and G^ are the same as in equations 
(IV-15) and (IV-13) . The eigenvalue equations are the same 
as equations (IV-16) and (IV-14) except that in (IV-16) h c 
is replaced with h£. Similarly, the slopes Aj^ and A^ 
are given by equations (IV-11) and (IV-12) after replacing 
h c with h£. 


Case C. 

In this problem, as in case A, the initial tempera- 
ture, t^, is uniform throughout the system. However in 
this case there are also heat transfer resistances at the 
outer boundaries characterized by the conductance coeffi- 
cients h, and h„ . The interface contact conductance co- 

efficient, h , is still present, thus giving gradient type 
c 

boundary conditions at all three physical boundaries. It 
should be noted that even though these solutions were devel- 
oped for the purpose of comparison with experimental results 
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for contact conductances at all three boundaries, they are 
equally applicable to the situation of convective fluid heat 
transfer at the outer boundaries with contact resistance in 
the interface. The forms of the boundary condition equa- 
tions, for the two phenomena are identical, and the only 
difference involved is the name given to h^ and h^. 

At time 0 — 0 , the x = 0 boundary is exposed to a 
medium at temperature t through the conductance coeffi- 
cient h^, while the x = L boundary is kept in contact with 
a medium at temperature t^ through the conductance coeffi- 
cient h^ • The contact conductance coefficient is kept con- 
stant at h^. In symbolic form these conditions are. 


T, (x,o) = O 

O * K £ Ol. 

(IV-35) 

T x C*,o) = o 

L 

(IV-36) 


x-<s> , e> o 

(IV-37) 


X-i , e >o 

( IV-38) 


x-cl , e>o 

( IV-39) 

Tj 

x-a , 0>c 

( IV-40) 


Equation (IV -8) is also applicable here and it is ob 
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vious that the same relation between the eigenvalues, i.e., 
equation (IV-14) , again occurs. 

As in case A the solutions may be written as a sum of 
steady-state and transient portions. The steady-state so- 
lutions are well known and similar to those of case A. As 
before, the transient solutions resulting from the separa- 
tion of variables in equation (IV-1) are written in a form 
which reduces the work required in obtaining the series co- 
efficients. 

Region (I) 

T|= K,00 + 21 t ***** ( iv -41) 

Region (II) 

— S & 

T x - K^cx) +■ ^ C " 1 


Here the steady state solutions are given by, 

K <*> = /- s (j£ 


(IV-42) 


wherein 


5 = 


i. + k. + + \ ^ 

Application of boundary condition (IV-37) gives 


(IV-43) 



(IV -44) 


It is noted that as 


becomes large. 


c n approaches zero 


84 


and thus agrees in the limit of an indefinitely large h 1 
with case A. 

Next equation (IV-38) is applied and results in 

/'■' jr f Si, LI 

U,T ’ - -I J (IV - 45 > 

where 

Nl* = (IV -46) 

As above, it can be seen that as h^ becomes very 

large , N approaches zero and the relation between G 
“ n 

and F n reduces to that of case A, as shown in equation 
( IV-13) . 


Application of the first interface boundary condition, 
equation (IV-40) , together with equations (IV-45) and 


(IV-14) gives 




AJ>, 4^. b — ^<>4 b 



(IV-47) 


Examination of equation (IV-47) shows that it also 
agrees with the limiting case. 

The eigenvalue equation results from the application 
of the second interface boundary condition. However, if the 
resulting equation is written out so that it involves only 
and the physical constants it is extremely involved. 
Since the results could never be used to any appreciable ex- 


tent without the aid of a digital computer it is simpler to 
leave the equation in terms of if and the grouped con- 
stants. This result is 



L 

*. 


C-n-UnY^ - 1J l /. 


(IV-48) 


Inspection of equation (IV-48) shows that it, too, reduces 
in the limit to the corresponding equation for case A. Only 


the series coefficients are needed to complete the so— 

lution. The method of obtaining the B^ is exactly as be- 
fore and the orthogonality factors, or discontinuous weigh- 
ing functions, are the same. The resulting equation for B n 
is slightly more complicated than before since two of the 
integrals, see equation (IV-19) , contain an additional term 
due to the fact that the C n is not zero. Several hand- 
written pages are required for the integration and simpli- 
fication, thus only the result is given here. 

Region (I) 



Region (II) 


( IV-49) 
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where E n is given by 


E h » ("Ct, 2 + +.2 -M*v a 4i* 
&)fe[ $21*: ~ Wuft lY/ , i c 

^ g„ 6 - A/* bj 1 ( N * H K 6 + 

S^b Cod fi^lo + 2 N ri -44vi. X S- h bJ 

This completes the solution for case C. 


( IV-50) 


Case D. 

In this problem the initial temperature distribution 

is linear, i.e., steady state heat flow. At time 9=0 the 

conductance coefficients are changed from h, , h and h 

1 c 2 

to new constant values h ' , h' and hi, while the outer 

1 2 

boundary medium temperatures, t Q and t^, are kept con- 
stant . 

The solution of this problem is identical to that of 
case C except for the numerator of the series coefficients 
B n . For this reason the steps of the solution will not be 
repeated here. The resulting temperature distributions are. 
Region I 

T«,«- l-5ft »■$)-£ e***'* 

71 = 1 ( IV-51) 

Region II 

T t w *(%<#)-£ 

2 risi ■*» 
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The eigenvalues ^ ^ are determined from equation 

(IV-48) by replacing h , h , h with h' h' , hi. For 

J- c 2 1 C 2 

the other terms equations (IV-44) through (IV-47) and 

(IV-50) may be used with the appropriate values of Y . 

w n 

The equation for S' is just like equation (IV-43) except 

that h' , h', h' should be substituted for h, , h , . 

x z 1 c 2 

The numerator term is given by 


H„=-2S 




( IV-52) 


It should be noted here that a solution for the prob- 
lem of case A has been previously found by Seide [192] . 
However, since Tittle's generalized orthogonality work [222] 
was not available, the solution of [192] was obtained by 
taking separate origins for the space coordinates of the two 
slabs and by judiciously choosing certain muliplying fac- 
tors. The validity of Seide 's solution is not doubted, but 
his procedure is not as straightforward as that used here, 
and the resulting solution is slightly more complicated. 

Thus only the form of the present solution for case A is 

since a prior solution of the problem exists. The 
solutions of cases B through D, for the more complex bound- 
conditions, are believed to be original contributions 
since no solutions of these problems were found in the lit- 
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erature. Actually only the solutions of cases C and D need 
be considered, since, as previously mentioned, cases A and 
B are special cases of C and D, respectively. 


Numerical Solutions 

To provide a check against possible errors in the ana- 
lytical solutions, a finite difference numerical solution 
technique was used to generate solutions for the same prob- 
lems. The nodal equations necessary for the numerical solu- 
tion may be derived either by applying finite differences to 
® equation and the boundary conditions, or by writing 
a simple heat balance for each node used in the solution. 
Both of these methods are well known and will not be repeat- 
ed here. A sketch representative of the one dimensional 
system for heat transfer purposes is shown in Figure 8. The 
following are the nodal equations which were used for the 
digital computer solutions. 

Region I 

For the x = 0 boundary node: 




(IV-53) 


For the interior nodes: 

T|« ~ + £. l (iv-54) 


L 
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-Finite Difference 
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For the contact boundary node: 

+[|-3M,-^7Ti£o ( IV-55) 


Region II 

For the contact boundary node: 
0 +-A® _ . . -T-® . „ Mi-e 


T,to**= 2 M Jatij+i +[' -2 M x -2 Tj“ , (IV-56) 


For the interior nodes: 

Uhi) = ^ t [ T i«-i)+1I(.„(]+[l-2MJX<„ pc-ni-Ht,} dv-57) 


For the x = L boundary node: 

+2M.T, v 0 +[ i-^;-2m,] t,«; 4 , 


(IV-58) 


where 


M, 


_ o( x A& 


N, = 


(AX) Z ' 

1 t' &X • 

. -v 6 - ■ 

o 11 , 4 * 


M„ = 


= 


N_ = 


(AX) 1 

L 

t c ** 




( IV-59) 


For those cases where there is no contact resistance at the 
outer boundaries, equations (IV-53) and (IV-54) are not re- 
quired. The equations for the outer boundary nodes would 
then be the same as for the interior nodes in the respective 


regions. 

In order to insure stable solutions it is necessary 
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that the coefficients of the last temperature in each of the 
equations (IV-53) through (IV-58) be positive. It was found 
that length increments of 0.1 inch and time increments of 
0.01 second achieved this stability in every case studied. 

Comparison of The Numerical and Analytical Solution 

The above nodal equations and the preceding analyti- 
cal solutions were programmed for use on a digital computer. 
For a number of cases, which represented the extremes in 
terms of thermal properties and hot and cold region lengths, 
both solutions were used for identical problems. The com- 
puter results were obtained in the normalized form indicated 
in the analytical solutions, i.e., 

~t (*,*) - ~ti 

to ~ ti 

In this form the "temperature" (actually it is a reduced 
temperature ratio) at any location, x, and any time, 0, is 
a number less than 1 and greater than zero. For the runs 
of the identical problems the results of both solutions were 
printed out for all the nodal locations (i.e., every 0.1 in.) 
and for times from 1.0 second to the steady state time. For 
reasons explained below the "steady state time" was taken to 
be the time at which the temperature drop across the slowest 
reacting portion of the system reached 99.0% of its steady 
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state value. These results were printed out by the computer 
to five decimal places. The results of the analytical and 
numerical solutions were always identical to the number of 
significant figures stated. Thus any difference between the 
results of the two solutions was less than one part in 10^. 
Since this accuracy was quite sufficient, no further checks 
were made to determine more exactly what the difference may 
have been. It should be emphasized that the solutions were 
checked in this way for the complete range of variables used 
in study. 

After the above-mentioned accuracy checks were made 
only the analytical solutions were used to generate the data 
for the study. The reason for this is that the analytical 
method requires less computer time. The nature of the pres- 
ent study was such that the desired output from the solu- 
tions was the "time to approach steady state" as defined 
above. With the numerical solution of this time dependent 
problem it is necessary to iterate the solution from time 
zero up to the steady state time for all the nodal points. 

In the present study, in which the desired times ranged from 
approximately 15 to several hundred seconds, such a process 
requires many iterations. On the other hand, the analytical 
solutions are "point" solutions. This means that it is pos- 
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sible to calculate the solution at any location and at any 
time. Advantage is taken of this fact on the computer by: 

(1) calculating the temperatures only at the desired loca- 
tions, i.e., at the boundaries and (2) by taking large time 
steps at first and then smaller steps to determine the 
steady state time. In the present study the savings in com- 
puter time was considerable. 

It should be noted, however, that the numerical solu- 
tions served the very useful purpose of providing a check on 
the analytical solutions. Also, if it were desired to vary 
the thermal properties or boundary conditions continuously, 
then the limitations of the analytical approach would neces- 
sitate the use of a numerical solution. 

Results of Computer Study 

The original purpose of the computer study was to fur- 
nish information for the design of experiments. These 
gedanken experiments provided data which aided in the selec- 
tion of metals to be tested, and which helped in the plan- 
ning of experiments with regards to test duration. It is 
for this reason that the 99 % criterion was chosen. It 
should be pointed out here that all the results presented 
below are for a situation corresponding to analytical solu- 
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tion cases A and B above. That is, for no contact resist- 
ance at the outer boundaries. At the time this work was 
done it was thought that the boundary conditions of the ex- 
perimental work could be made to approach this situation. 
However, it later proved impossible to achieve these condi- 
tions experimentally. Thus for the comparisons between the- 
ory and experiment which are presented in section VI the an- 
alytical solutions of cases C and D were used. 

During the course of the early computer study, made 
for experiment planning, interest grew in two ideas brought 
out by the study. The first of these was the possibility of 
making a general correlation on the time to reach steady 
state. Such a correlation would be of some practical impor- 
tance . The second was a phenomenon which was found to occur 
in certain instances and which eventually came to be called 
the "overshoot phenomenon." 

Because of time limitations most of the effort in this 
study was concentrated on the time to reach steady state 
correlation. The results of this correlation study are 
shown in Figures 9 through 13. All those figures are plots 
of a dimensionless time 


( IV-60) 
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versus a "thermal time-constant ratio" 

c~ cK.ywi ( iv-6i) 

for various values of an "inverse Biot number" 

./* m XZ (IV - 62) 

The correlations show, as reflected by the above quan- 
tities, the effects of geometry, material and contact con- 
ductances. 

Figure 9 represents all cases where the materials in 
both regions are the same. In the data used to plot the 
curves of Figure 9 the total length (a + b) ranged from 0.5 
to 8.0 inches, the length ratio (^/a) ranged from approxi- 
mately 0.1 to 10, and the metals ranged from stainless steel 
to pure aluminum. The latter gives a range of thermal con- 
ductivity of 10 to 117 Btu/hr . -ft . -°F , and a range of ther- 

2 

mal diffusivity of 0.15 to 3.33 ft. /hr. Contact conduct- 
ances were in the range 25 to 4000 Btu/hr . -ft . -°F . For 
these cases the ^ parameter is simply the length ratio 
( / a ) • The overall significance of this figure is that it 
makes it possible to estimate how long it would take a com- 
posite system to reach steady state after being subjected to 
a sudden temperature change on one side. It should be noted 
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that this time is independent of the overall temperature 
difference. With experimental data Figure 9 could also be 
used to determine contact conductance by measuring the time 
to reach steady state, for a system of known thickness and 
materials. These curves also show, through the^ -parameter, 
under what conditions the contact conductance would signif- 
icantly affect the response of the system. 

As an example of the use of the curves, consider a 
system consisting of one inch of aluminum in contact with 
two inches of aluminum through a contact conductance of 280 
Btu/hr . -ft? -°F . This gives X =0.5 and /3 = 5. From Fig- 
ure 9, (S) - 11.6 which gives 6 = 69.6 seconds. The same sys- 
tem, with a contact conductance of 2800, would have X - 0.5 
and = 0.5, which would give Q = 40.2 seconds. 

An interesting and somewhat unexpected phenomenon is 
demonstrated by the curves of Figure 9. Consider a system 
of the same material with the same contact conductance (280 

A 

Btu/hr. -ft.-°F) and same total length as the above example, 
but with equal thicknesses, a = b = 1.5 inches. Here X = 1.0, 
/& ~ 3.35, and the curves give (9 - 4.8. This results in a 
steady state time Q = 40.5 seconds. In other words the sym- 
metrical system reaches its steady state sooner. The impor- 
tance of the symmetry is further demonstrated by looking at 
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still another system which has the same materials and con- 
tact conductance as the above example, but which has equal 
thicknesses of a = b = 2 inches. The steady state time for 
this system is 69.0 seconds. It is seen that although this 
system has 25% more material than the 1 inch — 2 inch sys- 
tem its steady state time is approximately the same. It is 
this symmetry aspect which produces the minimum in the 
curves of Figure 9. These minima show that if T>1 the di- 
mensionless time is controlled by region 2 (the cooler side) 
and that region 1 controls when 'CKl. While affects the 
position of the curves, the minimum is almost independent of 
• 

For the cases in which the materials in the two re- 
gions are not the same it becomes necessary to present the 
information on two separate plots for each material combina- 
tion. For example, in Figure 10 an aluminum-tin system 
(aluminum in region 1 and tin in region 2 ) is presented, 
whereas in Figure 11 the materials are reversed. Figures 12 
and 13 show the correlations for a stainless steel-tin sys- 
tem. While it would be desirable to present all this infor- 
mation in one or two plots, it should not be surprising in 
view of the complexity of the system that this may not be 
possible. It may be possible to accomplish this through 
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finding another parameter and/or combining two or more para — 
meters. A possibility for an additional parameter, as indi- 
cated by the analytical solutions, is the ratio of the ther- 
mal conductivities. However, time did not permit further 
effort in this study. Prom the curves of Figures 10 and 11 
it can be seen that the aluminum tin system demonstrates the 
same behavior and effects as were discussed regarding Figure 
9. However, for these dissimilar metals the minima do show 
more of shift away from % - 1 asyS varies. The curves for 
the stainless steel-tin system, i.e., Figures 12 and 13 show 
a large change in the minima locations when the materials 
are reversed. In Figure 12, where the better conductor is 
in region 1, the minima show a shift with , but still oc- 
cur in the neighborhood of 'X = 1. When the materials are 
reversed the minima locations shift to the vicinity of 'X = 
0 . 1 . 

The overshoot phenomenon, mentioned above, represents 
the cases in which the temperature drop across the contact 
exceeds, or "overshoots," its steady state value during the 
transient portion before the steady state is reached. This 
phenomenon is illustrated in Figure 14 for a system consist- 
ing of aluminum and stainless steel with equal thicknesses 
of 1 inch for contact conductances of 200 and 1000 Btu/hr.- 
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Fig . 14. — Contact Temperature Drop Overshoot 
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A 

Figure 14 is a plot of the ratio of the contact 
temperature drop to the steady state temperature drop, as a 
function of time. Both curves approach a value of unity as 
the system approaches steady state. For this particular 
case the overshoot amounts to approximately 290% for the 
higher contact conductance value and 190% for the lower val- 
ue. It was found that for the cases studied the overshoot 
occurs for'£'>l, regardless of whether the same or different 
materials make up the system. Although no general correla- 
tion of this interesting phenomenon was attempted, it was 
realized that it is clearly a characteristic for the tran- 
sient behavior and was used in the comparison of the experi- 
mental and theoretical work in section VI. 

All of the above work is for the situation correspond- 
ing to case A, i.e., an initially uniform temperature and a 
sudden temperature rise at x = 0. Figure 15 illustrates 
the type of transient behavior that occurs for the situation 
of case B. In Figure 15 the temperature distribution is 
plotted for several values of time for an equal thickness 
(a — b = 1 inch) aluminum system. The system was initially 
experiencing a steady state heat flux with a contact conduct- 
ance of 1000 Btu/hr . -ft?-°F . The contact conductance was 
suddenly lowered to 300 Btu/hr . -ft? -°F. Figure 15 shows the 


(e 'x)jl 
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Fig. 15. — Effect of Suddenly Changing Contact Conductance 
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resulting temperature profiles for 1, 2 and 5 seconds and 
both steady states. Although a large number of such cases 
were studied for the purpose of designing laboratory experi- 
ments, no general correlation was attempted. 

A glance at the solutions for cases C and D above 
shows that they represent very complex phenomena. They are 
quite useful solutions and can be applied to any particular 
combination of materials for which information is desired. 
Indeed, they are used in section VI for the comparison of 
the experimental results. However, anything approaching a 
general correlation will require a great deal of study and 
considerable computer time. 

Summary of Theoretical Study 

Analytical solutions for the time dependent tempera- 
ture distributions for some problems of one -dimensional two- 
region systems separated by an interface with contact re- 
sistance have been derived. These solutions have been 
checked on a computer by a fine-network numerical technique 
and found to be accurate. The results of a computer study 
made with these solutions was presented in the form of di- 
mensionless correlations for the practical quantity, time to 
approach steady state. An illustration of the contact over- 
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shoot, which is recognized as a characteristic of transient 
response, was also given. 


V. EXPERIMENTAL PROGRAM 


Objective 

The objective of the experimental program was to de- 
sign and operate an apparatus which would provide test data 
on the transient response of one -dimensional, composite met- 
al systems with contact resistance when subjected to thermal 
transients. Primarily, this test data was desired for the 
purpose of checking the applicability of the theoretical so- 
lutions of section IV to real composite systems. To this 
end an attempt was made to create experimental conditions 
which approached, as closely as possible, the theoretical 
boundary conditions. 

Types of Experiments 

Basically the experimental program consisted of mea- 
suring the temperature distribution as a function of time in 
two metallic cylinders in contact while they were undergoing 
thermal transients. Figure 16 is a sectional drawing of the 
test section portion of the experimental apparatus. The 
other primary quantities measured were the axial force press- 
ing the two cylinders together and the temperatures of the 
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heat source and sink blocks. 

The thermal transients were produced by disturbing the 
one -dimensional heat transfer through the test specimens. 
Three test variables were employed in producing these dis- 
turbances: 1) the source block temperature; 2) the axial 

force on the contact surfaces; and 3) the environmental 
air pressure of the heat transfer system. Three different 
types of disturbances were studied. These consisted of 
holding two of the above test variables constant and varying 
the third rapidly from one fixed value to another, and then 
holding it constant. 

The basic features of each type of test using the 
above test parameters are described below. The various por- 
tions of the apparatus which are referred to in these des- 
criptions are shown in Figure 16. The equipment and exact 
procedures are described in detail later in this section. 

All of the different types of test were designed so 
that they could be run sequentially. Therefore, each type 
of test is designated by the word "phase" and a number which 
denotes its place in the standard test sequence. 

Phase 1. This test corresponds to the theoretical 
case C of section IV . The test began in a vacuum with a 
uniform temperature, in both specimens, equal to the sink 


112 


temperature. Than the sink block and samples were 
raised by the loading system so that the upper specimen was 
held against the hot source block with a constant force. 

Data were recorded while the system responded and sought a 
steady state. 

Phase 2. After the completion of phase 1, phase 2 was 
started by applying a step change in the force holding the 
specimens in contact. This test corresponded to the theo- 
retical case D with a sudden increase in contact conductance. 

Phase 3. This test consisted of reversing phase 2. 
After phase 2 was completed the contact force was suddenly 
dropped back to its original value. This test corresponds 
to the theoretical case D with a sudden decrease in contact 
conductance . 

Phase 4. When phase 3 had reached steady state phase 
4 was started by suddenly letting air back into the test 
chamber. The correspondence with the theoretical was the 
same as phase 2 except that the increase in conductance was 
achieved by adding air to the contacts at constant contact 
force. 

Phases 5 and 6. These tests were repeats of phases 2 
and 3 except that the samples were in air at atmospheric 
pressure instead of a vacuum. 
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The above brief descriptions demonstrate the basic 
philosophy of the experimental program and the correspond- 
ences intended to meet the stated objectives. Figure 17 is 
a photograph which provides an overall view of the basic 
test apparatus. Details of the test samples, equipment, and 
procedures are presented below. 

Test Specimens 

Test samples were constructed of 2024-T351 aluminum, 
Armco Iron, and type 303 -MA stainless steel. All samples 
were made of one -inch diameter bar stock and had nominal 
lengths of either one or two inches. Each specimen con- 
tained five thermocouples nominally placed at the per-cent- 
of -length values shown in Figure 18. Actual measured dimen- 
sions of all test specimens were recorded and are presented 
in Appendix A, in which other pertinent sample data are also 
shown. The procedures used to prepare the test specimens 
are given next. 

The sample material was first placed in a lathe, very 
lightly turned and then polished to make the cylindrical 
surface as reflective as possible. This was done to reduce 
both the absorption of radiation from the source block and 
the emission of radiation from the samples. The samples 
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Fig. 17. — Photograph of Test Apparatus 
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were then cut to rough length in a draw-cut saw. These 
blocks were then put back into the lathe and both ends were 
turned to a smooth tool finish, leaving the total length 
about 0.010 inch longer than desired. The samples were tak- 
en to a commercial grinding shop where they were ground to 
a surface roughness of about 100 microinches (rms) . The 
samples were then mounted in a shaper chuck, and the ther- 
mocouple slots were cut with a specially prepared shaper 
tool according to the scheme shown in Figure 18. After this, 
the samples were placed on a lapping table to be ground. 

All grinding was done by hand, starting with a number 240 
9 r it and continuing with successively finer grits until the 
desired surface finish was obtained. The "outside" end of 
ever Y specimen (the end that would be placed against a source 
or sink block) was ground in the final stage with an alumina 
powder. This resulted in surfaces with roughnesses of 3 - 6 
microinches (rms) . 

Initially samples were prepared by machining the con- 
tact surfaces in a shaper to produce surfaces with parallel 
waviness of regular wavelength. However, test results on 
the samples prepared in this way showed extremely distorted 
temperature profiles. This was caused by the highly convex 
test surface produced by the shaping process. The deviation 
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was large enough in one case to be seen by the unaided eye 
when two such surfaces were held together. Thus the grind- 
ing process was used to prepare all samples for which data 
are reported here. 

All surface roughness measurements were made using a 
Brush Surface Measuring System Model MS-5000. Surface re- 
cords were made and kept on each sample tested. The average 
surface roughness of each contact surface is given in the 
sample data in Appendix A. 

An attempt was made to measure the microhardness of 
the contact surfaces. However the surface roughness made 
it impossible to identify the edges of the indentations 
clearly. Consequently, microhardness measurements were made 
on polished flats on small cylinders of the sample materials. 
These measurements were made on a Tukon Model LL microhard- 
ness tester. Diamond Pyramid Hardness Numbers were obtained 
on the three sample materials for indenter loads of 0.50, 

0.75 and 1.00 kilogram. These data are also given in appen- 
dix A. 

Preparation of the test specimens was completed by in- 
stalling the thermocouples. As noted previously the ther- 
mocouples were installed in shallow slots rather than dril- 
led holes (see Figure 18). The idea for this type of instair 


118 


lation originated at the National Bureau of Standards in ex- 
perimental thermal conductivity work by Watson and Robinson 
[252]. The reasons for using slots are threefold. First, 
it is difficult to establish good thermal contact between a 
thermocouple bead and the bottom of a drilled hole. Sec- 
ondly, a drilled hole produces a larger interruption in the 
sample cross section. Finally there is more uncertainty 
about the axial location of the thermocouple junction in a 
hole. The first two problems above can have disastrous ef- 
fects on transient temperature measurements. Early in the 
test program a one-piece sample (no interface contact) was 
constructed with a number of holes drilled to the centerline 
on one side and an equal number of slots cut on the other 
side at the same axial locations. All the holes and slots 

instrumented with thermocouples and a phase 1 type test 
was run. The temperature readings of the hole thermocouples 

showed large lags behind the slot thermocouples, as much as 
. o 

30 F about half way to steady state. Although this result 
was proof enough in favor of the slots, further comment is 
in order. Since the slots were cut with a width slightly 
less than the diameter of the thermocouple wire a very tight 
fit resulted from forcing the thermocouples into the slots. 
This provides intimate metal -to -metal contact and eliminates 
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the contact resistance problem. The slots were small and 
completely filled with thermocouple metal thus reducing the 
interruption of cross section and making the change in ther- 
mal capacity very small. Also, the press fit nature of the 
slot type installation makes the accuracy of the bead loca- 
tion equal to the accuracy of the slot location; whereas a 
hole must be drilled oversize and thus produces an addition- 
al possibility of error. The slot locations were measured 
just prior to thermocouple installation on the micrometer- 
equipped sample stage of the Tukon Hardness tester. In this 
way slot locations (center of the slots) could be measured 
to ± .002 inch. Measured slot locations of all test speci- 
mens are shown in Appendix A. 

In the transient type experiments of interest here the 
heat flow through the specimen cannot be measured directly. 
In order to calculate the contact conductance it is neces- 
sary to evaluate the flux from the thermal conductivity of 
the materials and the measured slopes of the temperature 
profiles. Since the conductance values thus depend directly 
on the thermal conductivity it was necessary to have accu- 
rate data for the thermal conductivity. To provide this in- 
formation a program was carried out by Mr. D. R. Williams in 
which the thermal conductivities of the materials used in 


L 
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this work were accurately measured. The details of the mea- 
surement method are presented in Mr. Williams' thesis [250] . 
The thermal conductivity data were fitted with linear least- 
squares curve fits over the temperature range of interest in 
this study. The values of thermal diffusivity for the spec- 
imen materials were estimated from data available in the 
open literature for similar alloys. Due to lack of infor- 
mation these values were assumed to be constant. The ther- 
mal properties are shown, along with the other sample data 
in Appendix A. 

Temperature Measurement 

All temperature measurements were made using Chromel- 
Constantan thermocouples. The thermocouples were all made 
from the same spool of wire purchased from the Thermo-Elec- 
tric Company. Number 26 gauge wires (0.0159 inch dia.) with 
a polyvinyl insulation were used to form all the measuring 
junctions. A smaller wire than this would be desirable from 
the standpoint of conduction error, however, this size was 
the smallest that could be butt-welded efficiently. Each 
test required 14 thermocouples in all: 2 each in the sink 

and source blocks and 5 each in the two test specimens. The 
same source and sink block thermocouples were used for all 
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tests. Their measuring junctions were formed by twisting 
and soldering, and they were installed by soldering into 
one — eigth— inch deep drilled holes in the copper blocks, as 
indicated in Figure 19. 

Each test specimen had its own set of five thermo- 
couples which were not removed after the initial installa- 
tions. The specimen thermocouple measuring junctions were 
formed by butt-welding the Chromel and Constantan ends to- 
gether with a small resistance welding unit. This process 
is definitely an art. After some experimentation the pro- 
per power settings and method of preparing the wire ends was 
found. The resulting thermocouples had a good, welded meas- 
uring junction and no bead, thus providing a uniform diam- 
eter for installation in the slots. Installation was also 
found to be an art, and, after practicing on some dummy sam- 
ples, the following procedure was found to be successful. 
With the sample held in a small vise, the thermocouple ref- 
erence junction was centered over the slot and tamped into 
the slot. For the tamping a single-edged razor blade whose 
sharp edge had been ground off was held against the wire and 
struck with a light hammer. After the wire was inserted to 
the bottom of the slot the upper edges of the slot were fold- 
ed over on the wire to insure the security of the wire. 


Fi- 
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nally, the wires were wrapped around the sample and a short 
piece of heat-shrinking plastic tubing was shrunk over the 
two wires. At the point where the two wires met at the 
sample surface (opposite the slot) a small dab of silicone 
rubber cement was used to prevent relative motion between 
the wires and the specimen (see Figure 18) . 

This last step in thermocouple installation brings up 
an important point. Since transient temperature measure- 
ments were to be made, it was necessary to minimize the 
amount of extra material touching the samples, because ad- 
ditional material would add to the heat capacity and affect 
the transient response. For this reason the samples could 
not be insulated as they usually are in steady-state work. 

The thermocouple circuits are shown schematically in 
Figure 19. Since each sample had its own set of thermo- 
couples, all thermocouples were connected with gold plated 
socket -and -pin connectors near the base plate of the vacuum 
system to the lead wires from the reference junctions. Care 
was taken in the location of the connectors to insure that 
they were at a uniform temperature to prevent any additional 
thermocouple effect due to the dissimilar metals. The lead 
wires to the 32°F reference junction were Number 14 gauge 
Chromel-Constantan wires to reduce the electrical resistance. 
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The 10 thermocouples had electrical connection via the metal 
samples. No further electrically common connections could 
be allowed in order to prevent the thermocouple wires from 
forming current loops and thus giving erroneous voltages. 
Consequently, each thermocouple had its own separate refer- 
ence junction. These junctions were made of the 14 gauge 
wire and electrically insulated by plastic tubing and sil- 
icone rubber. All junctions were kept together by means of 
a rubber band and held in the middle of a large insulated 
container which was packed, top-to-bottom, with crushed ice 
and water. 

The requirement of simultaneous, transient temperature 
measurement precludes the use of a potentiometer for taking 
test data. All data were recorded on a Honeywell Model 1508 
"Visicorder" at paper speeds of 0.1 or 0.2 inch/second de- 
pending on the samples being tested. Thermocouples used in 
this way have a small electrical current flowing at all 
times. Thermocouple reference tables are based on a bal- 
anced voltage reading and could not be used in this applica- 
tion. It was thus necessary to calibrate all the thermo- 
couples of each specimen set over the range of temperatures 
to be incurred in the test. In order to provide the flex- 
ibility required by the variation of temperature range for 
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each thermocouple and to provide impedance matching for the 
galvanometers it was necessary to put a pair of variable, 
series-parallel resistors in each thermocouple circuit. 

These variable resistors were all located in a small metal 
cabinet which was located, electrically, between the refer- 
ence junctions and the galvanometer inputs. By adjusting 
both resistors in each circuit the desired galvanometer de- 
flection and impedance matching could both be achieved. 

Galvanometer ranging was accomplished by putting the 
specimens in position in the test apparatus, establishing a 
steady state heat flow through the samples and adjusting the 
ranging resistors to obtain the desired deflections. After 
the ranging was set the thermocouples could be calibrated. 
Calibration was performed using the samples and instrumenta- 
tion in the same configuration as they would be in for a 
test. This was accomplished by using a portable calibration 
bath which could be set up adjacent to the test fixture. The 
arrangement of this insulated bath is shown in Figure 20. 
With all thermocouple circuitry in the configuration it was 
to be used in during the tests, the samples were placed in 
the metal cans and covered with aluminum filings and insula- 
tion. The oil bath was constantly stirred by a paddle wheel. 
Bulk changes in the oil temperature were produced by flowing 
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steam through the lines for heating, and water through the 
lines for cooling (see Figure 20). When the desired temper- 
ature of a given calibration point was neared the oil bath 
temperature was controlled automatically by a Honeywell 
"Thermoniter" unit and two 150 watt electrical heaters. A 
set of calibrated ASTM standard thermometers, used to pro- 
vide temperature measurements over the full calibration 
range, were employed as indicated in Figure 20. To provide 
a check on the actual sample temperatures, and to provide a 
means of determining when the specimens had reached equilib- 
rium, the center thermocouples of each specimen were moni- 
tored on a potentiometer. This was accomplished by discon- 
necting the leads (of these two thermocouples) to the rang- 
ing unit and connecting the potentiometer leads in their 
place. Thus connected, the ranging resistors and galvanom- 
eters were completely removed from the circuits. The tem- 
P era tures of these two thermocouples were determined by con- 
verting the balanced voltage reading of the potentiometer 

a set of NBS conversion tables. The two readings were 
always within about 0.1°F of each other and usually were 
within the same tolerance of the thermometer reading. After 
the temperature readings had stabilized, the two center 
thermocouples were put back into the recording circuit and a 
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record was made of all the galvanometer deflections. This 
process was repeated until all the desired calibration 
points had been taken. 

Since there was an electrical current flow in this 
application the temperature versus deflection curves were 
slightly non-linear. Two typical calibration curves are 
shown in Figure 21. All deflections, for calibrations and 
test runs, were read from the oscillograph records with a 
scale ruled to 0.01 inch. The calibration readings were 
punched on cards which were processed by a computer program 
bo fit a least —squares second order polynomial through the 
calibration data of each thermocouple. These polynomials 
were used by the test data reduction program to calculate 
temperatures from the test data deflection readings. The 
solid lines in Figure 21 are the curve fits for the calibra- 
tion data shown there. Each set of test specimens was cali- 
brated prior to the running of its respective test series. 

Source and Sink Blocks 

During the experiments the test specimens were held 
between the source and sink blocks as shown in Figures 16 
and 17. These blocks were constructed of OFHC copper and 
were hollow with inlet and outlet ports. A detail drawing 
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Oscillograph Trace Deflection - inches 
Fig. 21. — Typical Thermocouple Calibrations 
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of the blocks is given in Appendix C. The sink block tem- 
perature was held constant during the experiments by flow- 
ing room temperature water from a large tank. Steam from 
the building heating system was used to maintain the source 
temperature constant; the temperature was adjustable over 
the range of 290-310°F by means of controlling the pressure 
using a throttling valve. During the experiments source 
temperatures of 300-304°F were used. The source and sink 
temperatures were measured by Chromel-Constantan thermocou- 
ples as indicated previously. These thermocouples were from 
the same spool of wire as the specimen thermocouples. Be- 
cause good accuracy found for this wire, as indicated above, 
the temperature -deflection calibrations of these thermocou- 
ples were made using the NBS table conversions. The proce- 
dure used was again that of disconnecting temporarily to use 
the potentiometer. The steam and water temperatures were 
varied over the small ranges of interest to provide these 
calibration points. Second order curve fits for these cali- 
brations were made and used in the manner described above. 

Test Fixture and Loading System 

The basic test fixture, shown schematically in Figure 
16, consisted of two 3/4-inch thick plates connected by 
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three 3/4-inch diameter rods. This basic framework was made 
of type 304 stainless steel. Detail drawings of all the 
test fixture components are shown in Appendix C. The two 
parallel plates were used to mount the sample holding and 
loading components, and were adjustable by means of the 
three threaded connecting rods to accommodate different to- 
tal specimen lengths. The source block was held by an 
adapter and thermally insulated to reduce heat loss to the 
test fixture. As shown in Figure 16, this adapter trans- 
mitted the applied force through a spherical joint to the 
upper plate, to which it was connected by three small heli- 
cal springs. This arrangement was used to insure that any 
small misalignment in the test column would not cause the 
contact surfaces to be loaded unevenly. A second, similar 
adapter was used to hold the sink block. This adapter rest- 
ed on the load cell (Lockheed Electronics Model WR75-025 
Load Washer) which in turn rested on the loading rod of the 
hydraulic cylinder. This test fixture assembly rested on 
the base plate of the vacuum chamber as shown in Figure 17. 
^ii instrumentation and fluid lines were passed through the 
base plate with standard vacuum type feed-throughs. The 
vacuum system was a Consolidated Vacuum Corporation Model 
CV-108 vacuum unit. 
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The force control system is basically a closed-loop 
electro-hydraulic servomechanism. Figure 22 is a schematic 
drawing of the main components of the system. A double end- 
ed hydraulic cylinder applies the force to the test speci- 
mens. The cylinder is supplied with a constant pressure, 
constant flow stream of hydraulic oil by a pump and flow 
regulator. A proportional type servovalve controls the 
pressure differential across the piston of the hydraulic 
cylinder and thus controls the force applied to the test 
samples. The cylinder and servovalve combination is a Moog 
Servocontrols , Inc. Model 1725G servoactuator . This servo- 
actuator was originally designed as an engine positioner for 
the Titan I missile and as such was made to operate using 
piston position for feedback. This system was modified to 
allow the use of applied force for feedback. 

As shown in Figures 16 and 19, the applied force is 
transmitted through the force cell. The force cell is a 
strain gage device. It is used as one leg of a resistance 
bridge so that the unbalanced bridge voltage is proportional 
to the applied force. This voltage is used as a feedback 
signal as indicated in Figure 22. By means of a potentio- 
meter an input signal is supplied to a summing junction at 
a polarity opposite to the feedback signal. 


The difference 
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Pig. 22. — Force Control System Schematic 
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between these two voltages is the error signal which is am- 
plified and used to drive the servovalve. Thus when the 
feedback is equal to the input the servovalve holds the de- 
sired force on the system. The input control consists of 
two potentiometers connected by a two position switch. This 
arrangement allowed the step change in input force to be ac- 
complished. Tests made on the system showed that the re- 
sponse time to the step change was of the order of millisec- 
onds, thus, for practical purposes the load changes were in- 
stantaneous . 

In addition to providing a feedback signal for the 
controller the force cell signal was also amplified and rec- 
tified for use in recording the force on the samples during 
the experiments. A complete set of electrical circuit dia- 
grams for the force control system are given in Appendix D. 
The force trace on the oscillograph was calibrated by plac- 
ing another load cell between the source and sink blocks and 
comparing the trace deflection to the output of this second 
load cell, which had been calibrated on a compressive test- 
ing machine. 

Experimental Procedure* 


The following is a brief outline of the procedures 
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used in operating the test apparatus to obtain the data for 
the various phases. 

After the thermocouple circuits had been calibrated 
the samples were removed from the calibration bath and 
cleaned. Each end of both specimens was cleaned with ace- 
tone and absolute ethyl alcohol. The outside ends of both 
samples were coated with a thin film of Dow Corning silicone 
grease to reduce the contact resistance between the speci- 
mens and the source and sink block surfaces. Then the spec- 
imens were placed in the test apparatus. 

Phase 1. The samples were pressed together between 
the source and sink blocks with a pressure of about 10 psi. 
By a valving arrangement cooling water was circulated 
through both blocks for about 10 minutes to bring the sam- 
ples to a uniform temperature. Next, the samples and sink 
lowered by the force system so that the upper 
(hot side) specimen no longer made contact with the source 
block. The vacuum bell jar was then lowered and test cham- 
ber evacuation started. While the system was pumping down 
the cooling water was kept flowing through the sink block 
but was shut off from the source block. Normally it re- 
quired about 20 minutes to pump the test chamber down to the 
range of 1-5 microns hg. This pressure range was used be- 
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cause it was low enough to make the interstitial gas con- 
duction small, and because it was easily obtainable with a 
mechanical pump. When the desired test chamber pressure was 
reached the oscillograph was started and the steam to the 
source block was turned on. About 20 seconds was required 
to bring the source block up to the desired temperature. 
Finally, an input signal to the force controller raised the 
sink block and test specimens, bringing the upper specimen 
into contact with the hot source block. All data were con- 
tinuously recorded on the oscillograph record throughout all 
phases of the tests. 

Phase 2. When phase 1 was completed the contact pres- 
sure was raised suddenly to a new constant value by flipping 
the two position switch on the force controller. The con- 
trol system maintained a constant pressure until a new 
steady state was reached. 

Phase 3. The contact pressure was returned to its 
original (lower) value by flipping the control switch back 
to the first position. This constant pressure was main- 
tained until the steady state was reached. 

Phase 4. The vacuum pump was turned off. With the 
force control system maintaining the same constant pressure, 
air was allowed to enter the system by opening a hand valve. 
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This produced an increase in contact conductance and the 
system sought a new steady state. The test data show that 
this was effectively a very rapid change. That is, the tem- 
perature profiles show a response like that of the phase 2 
tests. This provides experimental verification of the the- 
oretical predictions of Aaron and Blum that the "thres- 

hold" ambient pressure, above which the contact conductance 
little affected by ambient pressure, is very low. 

Phases 5 and 6. These phases were accomplished in the 
same manner as phases 2 and 3, respectively. However, the 
changes in contact conductance produced by changing the con- 
tact pressure were different due to the presence of air. 

Data Reduction Procedure 

As mentioned above, all data were recorded continu- 
ously on an oscillograph. The trace deflections of each 
measured quantity were manually read from the oscillograph 
records with a scale ruled to 0.01 inch. The time inter- 
vals at which data were read varied with the specimens being 
tested. During the early portion of each experiment smaller 
time intervals were used because changes were more rapid. 

As the rate of change became slower the "data times" were 
more widely spaced. Typical data time intervals were as 
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follows: every 2 or 4 seconds for the first 20 seconds, 5 

or 10 seconds for the next 40 seconds, 20 or 40 seconds for 
the next 120 seconds, and then every 50 or 100 seconds de- 
pending on the total time. The oscillograph was equipped 
with an automatic timing system which placed lines on the 
oscillograms every second. This system, which used an RC 
circuit, was checked and found to be slightly fast. For 
convenience the data were processed using the timing line 
time and then the small correction factor was applied to all 
the final results. The deflection reading of all 15 traces 
were read for each data time and recorded in a log book. 
These data were then punched onto cards to be read into the 
data reduction computer program. For all the test data re- 
ported here this amounted to approximately 2300 data cards. 

The data reduction computer program was written in 
FORTRAN and used on the Control Data Corporation's Model 
1604 and 3400 computers. The following is a brief outline 
of the steps performed by the data reduction program. 

1* Read in the required specimen data, such as ther- 
mal conductivity, exact sample length and thermo- 
couple locations. 

Read in the thermocouple and force trace calibra- 
tion data, i.e., coefficients of the least-squares 
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curve fits. 

3. Read in the experimental data cards. 

4. Calculate the temperatures and contact pressure 
at each data time using the calibration results. 

5. Print out all the above information in tabular 
form. 

6. Perform transient analysis. This consisted of 
making least-squares curve fits through the tem- 
perature profiles for each data time and extra- 
polating these to the boundaries to calculate the 
contact surface temperatures. Then the contact 
heat flux was calculated from the local slope of 
the temperature profiles and the thermal conduc- 
tivity. Finally, the contact temperature drops 
were calculated from the extrapolated profiles, 
and the contact conductance was determined from 
these temperature drops and the heat flux. 

7. Print out the results of the transient analysis 
in tabular form. 

The above outline provides an indication of the large 
amount of work necessary to reduce the data of an experiment 
of this type. It is estimated that manual reduction of the 
data presented in this work would have required several 


140 


months. Some further remarks are in order regarding item 
6 of the above program. The least-squares curve fits for 
the temperature profiles used in the program were second or- 
der polynomials. This form was chosen after a study was 
made in which several functional forms of curve fit were 
tried on some of the experimental data. The second order 
polynomial consistently gave the smallest average deviation 
of the forms tried. Higher order polynomials and other 
functional forms were tried but were found either to give 
larger deviations or to possess bad extrapolation behavior. 
For example, with five temperatures one might fit a fourth 
order polynomial and obtain zero deviation, but the curve 
may rise or fall sharply outside the fitted range and pro- 
duce meaningless extrapolations. 

As discussed later in section VI, the test results 
were evaluated on a basis of time to approach steady state 
or time to reach maximum overshoot. These times were deter- 
mined from hand made plots of the transient analysis. This 
was necessary because the temperature drop across a specimen 
versus time is an S-shaped curve making it difficult to use 
a curve fit. By hand plotting, smooth curves could be drawn 
through the data points for the extrapolated temperature at 
each end of a specimen versus time. The time at which 


a par- 
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ticular temperature difference occurred was determined by 
using a drafting divider. 

The test data and results obtained using the above 
procedures are presented and discussed in section VI. A 
brief evaluation of the accuracy of the experimental results 
concludes this section. 


Experimental Accuracy 

A discussion of possible sources of error in the pri- 
mary measured quantities is presented first. The effect of 
the errors on the calculated quantities is then evaluated. 

The procedure used to calibrate the contact pressure 
trace was described above. The slope of the pressure trace 
(psx per inch of deflection) was approximately 30 psi/inch. 


With a reading accuracy on the scale used of ±. .01 inch 
the contact pressure accuracy would be dt 0.3 psi. To this 
should be added the accuracy of the calibration device 
which was approximately ± 0.5 pound and therefore about 
0.7 psi in terms of contact pressure. Thus the contact 
pressures should be good to at least 1 psi. 

For the temperature data, that part of the possible 
error due to reading accuracy depends on the range setting 
of the galvanometer trace. The thermocouple nearest the 
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source had the largest deflection slope (°F per inch of de- 
flection) since it had to measure the highest temperature. 
This slope was approximately 45°F/inch which would give a 
reading accuracy of ± 0.45°F. The thermocouple nearest the 
sink had a slope of about 30°F/inch, or a reading accuracy 
of ± 0.3°F. The average deviation of the calibration curve 
fits from the calibration data was about cb 0.3°F. The tem- 
perature profile curve fits showed average deviations from 
the data which varied with time from about 0.6 to 0.2°F, the 
former occurring at early data times and the latter as 
steady state was neared. Assuming an average value of 0.4°F 
for this deviation and assuming that the probable error is 
the square root of the sum of the squares, the probable er- 
ror in a calculated temperature is around ± 0.6°F. Thus a 
calculated temperature difference is probably only good to 
about rfc 1°F . 

The contact conductance was calculated by dividing the 
contact heat flux by the contact temperature drop. The heat 
flux was calculated as the product of the local temperature 
profile slope and the thermal conductivity. Actually, the 
average of this product on both sides of the interface con- 
tact was used in these calculations. Symbolically, 
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The thermal conductivity data is believed to be within 


(V-l) 
± 4 % 


^250] . It is practically impossible to evaluate the accu- 
racy of the temperature profile slopes during the transients 
and it can only be estimated in the steady state. However, 
only steady state values of conductance were used in the cor- 


relation of test results presented in section VI. In the 


steady state the temperature profiles are essentially linear. 
Assuming the worst situation, namely all the temperatures in 
one half of a specimen are high and those in the other half 
are low, would give an error of about 0.5°F/inch in the 
slope. The lowest slope measured was about 4°F/inch but it 
was usually around 20-30 °F/inch 

The percentage accuracy of the calculated contact con- 
ductance depends somewhat on the level of conductance, pri- 
marily because of the contact temperature difference (see 
equation V-l) . As an illustration consider the steady state 
value of phase 1 of run 405-2. This represents the lowest 
value of interface conductance calculated in these experi- 
ments. The contact conductance was calculated as 29 Btu/hr.- 
2 

^ ~ ^ * The contact A T was 142. 9°F and the slopes were 
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8 . 6°F/ inch on the hot side and 76.7°F/inch on the cold side. 
Thus the possible percent error in A T is about 0.6% and the 
possible percent error in the slopes were about 5.8% and 
0.6%. All errors combined would give a result of about ±10% 
for the contact conductance. For this same run the hot end 
conductance (conductance for the contact between the source 
block and the upper specimen) was calculated as 848 Btu/hr.- 
from a A T of 4.8°F. The possible error in AT is 
now about 20% giving a conductance accuracy of about 25%. 

The highest value of interface conductance calculated 
in the present data (2505 Btu/hr . -ft? -°F) was on phase 5 of 
run 304-8. The contact temperature drop was 21.0°F, the 
slopes were 28 and 60.5 °F/inch. The probable error in AT 
would be about 5% and for the slopes about 2% and 1%. This 
again results in about ±10% for the conductance. For this 
experiment the hot end conductance was calculated to be 5063 
Btu/hr . -ft . -°F . with a AT of 10.4°F, or, in other words, 
about ± 15% for the conductance. 

Based on this approach, it is believed that the con- 
ductance values are good to about 10% at best and probably 
average closer to 15%. It should be noted, however, that in 
some cases on the 405 series the contact temperature drop on 
(grease filled) end contacts was so low that the calcu — 
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lated conductances exceeded computer print format. In some 
of this series even negative values were calculated for the 
end conductances, which demonstrates the problem of A T ac- 
curacy. However, a check of the data shows that any calcu- 
lated end conductance below about 5000 Btu/hr . -ft?-°F is 
probably at least within ±. 20 %. 

Although these estimates of data accuracy certainly 
leave a lot to be desired one is hard-pressed to find many 
experimental contact conductance works that show a better 
accuracy. Some suggested changes for possible future work 
which might improve upon the overall data accuracy are of- 
fered in section VII. 


VI . EXPERIMENTAL RESULTS 


The experimental data obtained in the present work are 
presented in this section. All data were obtained employing 
the equipment and procedures described in the previous sec- 
tion. As used in this section, the word M data" refers to 
the printed output of the computer data reduction program. 
This data consisted of the following: (1) temperature dis- 

tributions in the two test samples and the source and sink 
block temperatures; and (2) the quantities which were cal- 
culated from the temperature distributions and thermal pro- 
perties, all at specified time intervals. The latter group 
includes the contact temperature drop, ^T , and the three 
contact conductances h c , h^ and h^ (see Figure 7). These 
data are presented below in four separate groupings accord- 
ing to the important point emphasized and the type of test. 

Contact Conductance as a Function of Time 

One of the first effects noticed in analyzing the test 
results was the apparent variation of the contact conduct- 
ance coefficient with time. A plot is shown in Figure 23 of 
the calculated contact conductance coefficient versus time 
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Fig. 23.-- Calculated Contact Conductance as a Function of Time 
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for two typical test runs. The sharp variations shown in 
Figure 23 were unexpected, and it was decided to attempt to 
rmine whether such variations were real or whether they 
were a result of the calculation procedure. For this per- 
pose several sets of data were generated with the theoret- 
ical solution for a constant contact conductance. These da- 
ta consisted of the theoretical temperatures at the loca- 
tions in the samples where the temperatures were measured in 
the experiment, at the time intervals used in the experi- 
ments. Next, these data were punched on cards and fed into 
the data reduction computer routine just as if they were ex- 
perimental test data. The results of two such check runs 
are shown in Figure 24. It can be seen that even though the 
data used were generated using a constant contact conduct- 
ance, the data reduction program calculates a conductance 
which varies with time in much the same manner as shown in 
Figure 23 for experimental data. As a further check on this 
effect, some of the same data was run through another compu- 
ter program written by Dr. James Beck of Michigan State Uni- 
versity [249]. Dr. Beck's program employs a different meth- 
od of calculating the conductance coefficient. However, the 
results were the same. That is, this program predicted the 
same kind of variation for both experimental and theoretical 
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24. — Check Runs With Theoretical Data 
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data. This evidence is not conclusive proof that the indi- 
cated variation in contact conductance is a result of the 
calculation procedure, but it does show that such a varia- 
tion is calculated when none exists. 

There is nothing in any of the theoretical work on 
contact conductance [55,58,85,143: to suggest a variation as 
strong as the one indicated in Figure 23. The only varia- 
tion that might be suggested by the theories of contact con- 
ductances is a slight variation due to the changing mean in- 
terface temperature. Figure 23 shows that after the initial 
sharp variation has settled out there is a gradual rise. 

This rise is of the order of that found by others studying 
the effects of mean interface temperature [24,27,47,140:. 

Of the reported works which have measured contact con- 
ductance by a transient means [26,46,60,127,190,231: only 
one [190: h as indicated a variation of the order indicated 
in Figure 23. The others have either not mentioned observ- 
ing any variation [60: , found it to be of the nature of scat- 
ter [26,46,127:, or found that it was a result of mean tem- 
perature increase [ 231 : • Schauer and Giedt [190: measured 
the contact conductance between two thin strips during tran- 
sient heating. Their method and results were discussed in 
section iv. They found a sharp increase in contact conduct- 


151 


ance which occurred over a period of approximately 100 mil- 
liseconds. The present writer feels that their results may 
be questionable. However, assuming that they are valid 
would suggest that any sharp variations in contact conduct- 
ance would have disappeared long before any time at which 
data were recorded in the present experiments. 

In view of all of the above evidence it is the writer's 
opinion that the contact conductance variation indicated in 
Figure 23 is not real but is a result of the calculation 
method. This can be explained as follows. During the early 
data times the contact heat flux is small and the resulting 
contact temperature drop, AT c , is small. This results in 
low accuracy for the calculated contact conductance because 
the error incurred in extrapolating the curve fits of the 
temperature profiles is not small compared to the calculated 
AT C . in fact, in some cases, the early-time conductance 
values were calculated to be negative as a result of the low 
accuracy in the contact temperature drop. Another source of 
error in the early conductance calculations is the tempera- 
ture profile slopes. Apparently the early-time temperature 
profiles vary too sharply near the contact plane to be close- 
ly fit by the curve fit used. As a part of the data reduc- 
tion routine the average deviation of the temperature pro- 
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files were also calculated at each data time. The devia- 
tion was always worse during the early times and improved 
for later times. For example, the data of the phase 1 por- 
tions of the 304 series showed deviations of around 0.6°F 
at the early times and decreased to around 0.1°F for the 
later times. The slope accuracy problem is made worse by 
the fact that the cold-side specimen temperature profile is 
slower in developing because of the contact resistance. 

Even after the calculated temperature drop across the con- 
tact has become sufficiently large (to prevent its being 
calculated to be negative) this second source of error may 
^e present. As the heat flow increases, the calcu- 
lated AT c becomes more accurate and the profiles become 
more amenable to the curve fitting. This results in better 
contact conductance values as evidenced by the above-men- 
tioned check runs of Figure 24. From an "information theo- 
ry point of view the above remarks might be summarized by 
saying that during the early times, when the heat flux is 
low, the five thermocouples in each specimen do not provide 
su ^^ c ^ en ^ information to calculate the contact conductance. 
This suggests that better values might be obtained by em- 
ploying more thermocouples in each specimen to provide a 
better knowledge of the actual temperature profiles. As 
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discussed in the previous section, it was sometimes possi- 
ble to use only four of the five thermocouple readings in a 
test specimen due to the distortions caused by the contact 
surfaces. in view of the above discussion this probably 
compounded the early time conductance calculation problem. 

In summary it is stated that the small gradual rise 
m calculated contact conductance at the later times, indi- 
cated in Figure 23, is believed to real. The sharper var- 
iations at the earlier times are believed to be only a re- 
sult of the calculation method. Thus, the evaluations and 
comparisons of the data which are presented below are based 
on the contact conductance values calculated when a steady 
state had been reached. 

Co mpariso n Wit h th e Results of Other Investigators 

As a means of providing a rough, overall check on the 
accuracy of the data obtained in this study, the values of 
contact conductance as a function of contact pressure were 
compared to results published by other investigators. These 
comparisons are shown in Figures 25 and 26. 

In Figure 25 the vacuum data of the 304 series are 
compared to some data obtained by Fried [922 and Clausing 
and Chao C58] . The specimens used in the 304 series were 
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Fig. 25. — Comparison With Data of References 58 and 92 
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made of 2024-T351 aluminum with surface roughnesses of 70 
and 75 microinches (rms) . The cited data of Fried were ob- 
tained for 2024-T4 aluminum for surface roughnesses of 45 
and 50 microinches (cla) , and that of Clausing and Chao were 
for 2024-T4 aluminum for 45 and 80 microinches (rms) . Based 
on the surface roughnesses the comparison indicates that the 
present data are reasonable. The absence of data on flat- 
ness deviation for the present data and Fried's data make a 
more valid comparison impossible. However, since the accu- 
racy of the present steady-state data, which was discussed 
in the previous section, is at least as good as that of the 
other investigators, and since considerable scatter is known 
to exist in all the literature Ql03,12^ , the above compar- 
ison is adequate to demonstrate that the present data are 
reasonable . 

In Figure 26 the in— air results of the 7A4 series are 
compared to the data of Fenech and Rohsenow [85] . The spec- 
imens of the 7A4 series were 2024-T351 aluminum and Armco 
Iron with surface roughnesses of 30 and 75 microinches (rms), 
respectively. The data of Fenech and Rohsenow are for spec- 
imens of aluminum (alloy unspecified) and Armco Iron with 
surface roughnesses stated to be 150 microinches maximum (no 
average given) . The comparison shown in Figure 26 is, again. 
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Fig. 26. — Comparison With Data of Reference 85 
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only offered as an indication of reasonableness. 

As indicated in the above discussion, it is very dif- 
ficult to make a better comparison of the test data because 
of the impossibility of matching even the most important 
parameters of materials, surface roughness and contact pres- 
sure range . 

Before proceeding to the main point of interest in the 
present work, there is one other observation which bears 
mentioning. The data of series 304, shown in Figure 25, il- 
lustrate quite clearly the hysteresis effect mentioned ear- 
lier. Going from a given contact pressure to a higher con- 
tact pressure and then returning to the original value again 
is seen to result in a higher value of contact conductance 
the same pressure. This result has been observed by oth- 
ers [85,92,233] but their tests were always run by continu- 
all y increasing and then continually decreasing the contact 
pressure. The present results show that even when the ini- 
contact pressure is well below the maximum value to 
which the surface has been subjected, the effect is still 
present. 

The vacuum data of the 7A4 series are also shown in 
Figure 26. The effect of "breaking in" the surface is dem- 
onstrated in these data. These data are for a set of spec- 
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imens whose surfaces had been exercised in previous tests, 
whereas the data shown in Figure 25 were for the first tests 
made on those surfaces. Comparing Figures 25 and 26 shows 
that the hysteresis effect was considerably reduced by hav- 
ing the surfaces well broken in. This affords an example of 
another variable in contact conductance - surface loading 
history. 

Comparison of Experimental Results with Theory 

The primary purpose of the present study was to deter- 
mine whether the one-dimensional theoretical solutions ob- 
tained in section IV could be used to predict, with suffi- 
cient accuracy, the transient behavior of a one— dimensional 
two }-ayej; with contact resistance in the interface. 

To this end, a series of tests were run in which the bound- 
ary conditions were made to closely approximate those of the 
theoretical solutions. The results of these experiments are 
compared to the corresponding theoretical cases below. 

Phase 1 Results 

As indicated in the discussion in section IV, the time 
to approach steady state is believed to be a parameter which 
is characteristic of transient behavior, as well a practical 
For this reason the time to approach steady state is 
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one of the parameters chosen as a basis for comparing the 
experimental results with theory. For use in those compar- 
isons the time to approach steady state was defined to be 
the time at whiph the temperature drop across the slowest 
reacting portion of the system was within one "time con- 
stant" of its steady state value, i.e., to within e - ^, thus 
the fraction is 1-e -1 ^ .632. It was necessary to go to a 
lower value than the 0.99 fraction used in the theoretical 
correlation to reduce the sensitivity to experimental error. 
The 99% criterion is more meaningful from a practical stand- 
point, but in this regime the rate of change of temperature 
drop with time is very low and small errors in temperature 
measurement would result in large time errors. 

The results of the comparisons of the Phase 1 data are 
presented as plots of time to approach steady state versus 
interface contact conductance. The actual experimental val- 
ues are tabulated in Appendix B. Those values shown in the 
plots reflect slight adjustments in the times made as fol- 
lows. For each run of a given series, i.e., set of spec- 
imens, the contact conductances on the ends of the speci- 
mens (the outer boundaries) were different due to the dif- 
ferent contact pressures on each run. It was desired to 
show the data in a form that could be compared to a single 



160 


theoretical curve for each series, rather than a point-for- 
point comparison which could not be meaningfully plotted. 
The test data were plotted on a "map" of theoretical data 
which consisted of a plot of time to approach steady versus 
interface contact conductance for various values of the end 
conductance. From this plot it was found that the experi- 
mental data exhibited about the same slopes as the theoret- 
ical for dependence on the end conductances. Thus the ex- 
perimental times to reach steady state were adjusted by ad- 
ding or subtracting a small amount. This amount is the dif- 
ference in time found from the theoretical curves, at the 
experimental value of interface conductance, in going from 
the end conductance measured in each experiment to a value 
which was the approximate average for the complete series. 

In this way experimental values of steady-state time versus 
interface conductance for the same end conductance could be 
determined for each test series and compared with theory. 

It is emphasized that, as can be seen in Appendix A, the ad- 
justments were usually small and should not affect the va- 
lidity of the comparison. This is indicative of the fact 

that when the end conductances, h and h are large com- 

1 2 

pared to the interface conductance, h c , they do not exert as 
large an influence on the time. This was the case in all 
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the present data. 

The resulting comparisons made in the above -described 
manner are shown in Figures 27 through 31 for the six sets 
of samples tested. In these figures the solid lines repre- 
sent the theoretical predictions made using the theoretical 
solution corresponding to Case C of section IV. A constant 
thermal conductivity calculated from the experimental data 
[|25CQ at the average specimen mid— point temperatures 
used in the theoretical solution. The experimental data 
points are plotted with symbols. It can be seen that the 
agreement between theory and experiment is better for some 
test series than for others. The worst comparison is for 
the 304 series, shown in Figure 27. The agreement for the 
904 series is somewhat better, as seen in Figure 28. Figure 
29 shows the comparison is still better for the 407 and 47A 
series, both on a straight difference and percentage differ- 
ence basis. The best agreement was found for the 7A4 and 
405 series, as shown in Figures 30 and 31. It is obvious 
that the agreement is better for the slower reacting systems. 
It is believed that this can be explained by the failure of 
an experimental boundary condition to match the theoretical 
boundary condition - specifically, the step rise in the 
source temperature, it was observed that during the experi- 
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Fig. 27. — Time to Approach Steady State for Series 304-Phase 1 
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Fig. 28. — Time to Approach Steady State-Series 904 Phase 1 
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Pig. 29. — Time to Approach Steady 
State-Series 407 and 47A Phase 1. 
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ments the temperature of the source block (copper block- 
steam system) would "dip” slightly when the upper specimen 
is brought into contact with it because of the finite capac- 
source block. The dip amounted to about 4 — 6°F 
at the maximum and required approximately 10 - 30 seconds to 
vanish completely, depending on the system. In other words 
the overall driving force ( total A T) was low for the early 
part of each experiment. As can be seen in Figures 27 - 31, 
and as would be expected, the error thus produced is largest 
for the systems with the lowest total resistance, which 
would also have the lowest time to reach steady state. This 
point is illustrated by comparing the Figures 30 and 28, the 
former is for a 2-inch Armco Iron specimen above a 2-inch 
aluminum specimen (series 7A4) , where as the latter is for 
two 2-inch aluminum samples (series 904). The Armco Iron 
sample being a poorer conductor than aluminum does not cause 
as much of a dip in the source temperature. The dip also 
does not last as long and this time represents a smaller 
portion of the time to reach steady state. Thus the error 
produced is much smaller. Similar arguments explain why the 
agreement was better for two 2-inch aluminum specimens. Fig- 
ure 28, than for a 1-inch above 2-inch aluminum set. Figure 
27; as well as why the Armco Iron-aluminum specimens. Figure 
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30, compared better than the reversed aluminum-Armco Iron 
specimens. Figure 29. For the 2-inch aluminum above 2-inch 
stainless steel series, Figure 31, the dip was about the 
same as for the 904 series, but the time represented a much 
smaller part of the steady state time, which resulted in 
better agreement. The result of this argument is that the 
agreement between theory and experiment is good for those 
experiments which more closely matched the theoretical 
boundary conditions. 

It should be emphasized here that Figure 29 contains 
the results of two test series, 407 and 47A. These series 
used the same specimens except that the surface of the Armco 
Iron specimen was re-ground to produce a different surface 
roughness (see Appendix A). Also, as indicated in Figure 
29, two of the runs for the 407 series were made in air. 

This provided two additional means of varying the interface 
contact conductance, h c , and thus furnished a further check 
on the validity of h c as a correlating parameter. No effect 
due to surface finish or the presence of air was noted. 

For another comparison between theory and experiment 
the overshoot in the interface contact temperature drop was 
used. As explained in section IV, under the condition that 
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T= s# >/ 


(VI-1) 


the temperature drop across the interface contact ( 4 tJ ex- 
ceeds or overshoots its steady state value for a while dur- 
ing the transient period. The quantity used for comparison 
in this phenomenon is the time required for the transient 
AT C to reach its maximum value. In the theoretical work it 
was found that this time depended on h c . The experimental 
values of this parameter are shown in Appendix B. Small ad- 
justments were also made to this parameter, in the manner 
discussed above, to provide a set of experimental data for a 
single value of the end conductance, h^ and , for each 
series. The results are plotted in Figures 32 through 34, 
in which the theoretical curves are drawn as solid lines and 
the experimental points are plotted as symbols. It can be 
seen that the agreement follows about the same pattern as in 
the time to approach steady state comparison. Note that the 
error is such (as in the previous case) that the experimental 
observations appear to be somewhat slower than the theoret- 
ical. That is, the experimental time of maximum overshoot 
occurs later than the theoretical time. This also seems to 
support the above argument in that a temporary dip in driv- 
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Fig. 32. — Time of Maximum Overshoot -Series 304 Phase 1 
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ing potential should have a delaying effect on the time to 
reach maximum overshoot. As in Figure 29, Figure 33 repre- 
sents tests with two surface roughness combinations and 
tests in air and vacuum. Again no effects due to these 
quantities were noted. 

As further evidence in the comparison it is pointed 
out that only four of the six series produced an overshoot 

ln AT c* For the tw ° series in which r < 1 (904 and 7A4) 
no overshoot was found to occur, which agrees with the the- 
oretical work. 

The phase 1 results may be summarized as follows. The 
poorer comparison of the 304 and 904 series is attributed to 
the source temperature discrepancy explained above. For the 
407-47A (same samples), 7A4 and 405 series the agreement be- 
tween experiment and theory is considered to be very good, 
especially in view of the experimental complexity. 

Results of Test Phases 2, 3,4, 5, and 6 

After each phase 1 experiment had reached its steady 
state a sequence of other tests were made with the samples 
still in place. The details of the test procedures used were 
discussed in section v. Each of these tests corresponds to 
the theoretical Case D of section IV. That is, the system 
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has a steady state heat flow established, and at time zero 
the interface conductance, h c , and the end conductances, 
and h 2 , are suddenly changed to new, constant values. 
Each succeeding test was started after the preceding phase 
had reached its steady state. For phase 2 the change was 
produced by suddenly increasing the contact pressure and 
holding it constant. Phase 3 consisted of suddenly decreas- 
ing the contact pressure back to the original value (the val- 
ue at the end of phase 1) . The change was produced in phase 
4 by suddenly letting air into the vacuum chamber. Phases 5 
and 6 were the same as phased 2 and 3 except that the con- 
tact pressure changes were made with the specimens in air 
instead of a vacuum. This sequence of tests provided exper- 
imental data for comparison in which the changes in conduct- 
ance were both increasing and decreasing, and which were pro- 
duced by two separate means . 

For this second type of test the parameter used for 
comparison of the experimental data and the theoretical so- 
lutions was the time to approach steady state; again, be- 
cause it is an obvious and practical characteristic of tran- 
sient behavior. However, in these cases the dependence of 
the time to approach steady state on the contact conduct- 
ances at the three broudaries is much more complicated than 
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for the phase 1 case. The theoretical solution in section 
IV (Case D) shows that transient portion — and thus the 
time to approach steady state — depends not only on the 
final values of h£, h^ and V (values after change) but 
also on the old values h c , ^ and h . It is apparent that 
such a dependence prohibits the results from being displayed 
in as simple a manner as was possible for the phase 1 data. 
Hence, it was decided, due to the complexity of the situa- 
tion, that the only practical comparison that could be made 
was a pomt-for-point comparison. That is, to compare the 
experimental value of the time to approach steady state with 
the value predicted by the theoretical solution for the ex- 
act conditions of each individual run. The theoretical val- 
ue for each run was found by putting the experimental values 

of the steady state contact conductances, h h' h h' 

c f c' n 2 

and h 2# for that run into the computer programmed solution 

and calculating the theoretical time to approach steady 
state . 

When the above scheme was first employed the same cri- 
terion for the time to approach steady state as was used in 
the phase 1 was used, namely, 1-i , or 0.632. It was 
found that the resulting comparisons were very poor. The 
reason for this can be explained as follows. First, the 
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time to approach steady state for these cases is much small- 
er than for the phase 1 type transients, as can be seen in 
the tabulated data. For, example, the theoretical times to 
approach steady state for the 904 series (two 2-inch alumi- 
num specimens) were all below 1.5 seconds using the 0.632 
criterion. Secondly, a problem arises in curve fitting the 
early time temperature profiles as it did in the phase 1 
case. The problem is that when the contact conductances 
are changed suddenly the temperature distribution is the vi- 
cinity of the contacts responds very quickly, thereby pro- 
ducing a sharp change in temperature profile near the bound- 
aries. Thus the curve fit, which cannot match this change, 
produces erroneous boundary temperatures when extrapolated. 
This means that if the time to approach steady state occurs 
quickly enough it is bound to be erroneously determined. To 
circumvent this problem it was decided to go to a fraction 
larger than 0.632 for the criterion. This allows the tem- 
perature profiles time to smooth out so that the extrapola- 
tion of the curve fits to the boundaries more accurately 
represents the true boundary temperatures. A fraction of 
0.800 was arbitrarily chosen. It was felt that anything 
larger than this might cause accuracy problems due to the 
small temperature -time slopes which occur as steady state is 
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approached. Since there is nothing to suggest that there is 
anything "holy" about any particular number this change of 

cr ^^ er ^ on should not have any effect on the outcome to the 
results. 

Therefore, the experimental and theoretical times to 
approach steady state were all recalculated on the basis of 
the 0.800 criterion for all the test runs. The complete 
test results are presented in Appendix B. Since a point - 
for-point comparison was made the results pertinent to the 
comparison for these data are presented in Tables 4 through 
6. In these tables the test runs and phases are identified 
and the values of the initial and final conductances are 
< ^ v ® n * The experimental and theoretical times to approach 
steady state and the difference between them are also given 
in the tables. The difference shown in each entry is the 
experimental value minus the theoretical value. As can be 
seen the comparisons show both plus and minus differences. 
This indicates that the differences are of a scatter nature, 
in contrast to the phase 1 results, which tended to be off 
in one direction. This observation adds further evidence in 
support of the argument given in the discussion of phase 1 
results regarding the source temperature dip; because, in 
these phase 2, etc., 


experiments the temperature of the 
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TABLE 4 . 

Comparison for Series 304 & 7A4 




Contact Conductances 
(Btu/hr-ft -°F) 

Steady 
St. Time 



1 Interface 


■IHb 

Cold 

End 

IE33I 

The. 

Run 


In' 1 


In’ 1 

mm 

In' 1 

ERiM 

Sec . 

Sec . 

304-5 

2 

638 

2319 

3523 

4900 

3402 

4500 

5.7 

5.44 


3 

2319 

1238 

4900 

4169 

4500 

4195 

8.7 

9.40 


4 

1238 

1975 

4169 

4366 

4195 

4305 

7.2 

7.08 

304-6 

2 

199 

1707 

1200 

4444 

1838 

4366 

9.5 

8.08 


3 

1707 

446 

4444 

1993 

4366 

2923 


13.74 


4 

446 

1092 

1993 

2345 

2923 

3738 


6.42 


5 

1092 

2441 

2345 

4220 

3738 

4606 

9.2 

10.60 

304-7 

2 

389 

1893 

1972 

4932 

3225 

4747 

mm 

9.30 


3 

1893 

907 

4932 

2818 

4747 

3825 


11.74 


4 

907 

1549 

2818 

3510 

3825 

4077 

9.7 

10.38 


5 

1549 

2637 

3510 

4745 

4077 

4917 

8.7 

7.08 

304-8 

2 

1111 

2987 

4571 

5925 

4904 

5679 

7.2 

6.76 


3 

2987 

1657 

5925 

4697 

5679 

5400 

9.2 

10.50 


4 

1657 

2505 

4697 

5063 

5400 

5431 

6.1 

9.24 


5 

2505 

3553 

5063 

5776 

5431 

6007 

4.4 

5.64 

7A4-3 

2 

102 

640 

952 

4870 

1063 

3408 

27.1 

21.08 


3 

640 

121 

4870 

2123 

3408 

818 

53.0 

60.68 


4 

121 

408 

2123 

2399 

818 

1077 

32.9 

37.09 


5 

408 

1351 

2399 

4022 

1077 

3716 

37.6 

37.37 

7A4-4 

2 

117 

868 

1184 

3357 

1127 

4129 

28.1 

24.77 


3 

868 

146 

3357 

1771 

4129 

993 

53.2 

56.98 


4 

146 

498 

1771 

1741 

993 

1297 

35.8 

37.91 


5 

498 

2032 

1747 

3763 

1297 

4505 

24.2 

23.20 

7A4-5 

2 

196 

686 

1457 

2417 

2998 

5101 

25.1 

22.32 


3 

686 

234 

2417 

1752 

5101 

3481 

38.8 

43.14 


4 

234 

808 

1752 

1725 

3481 

4055 

31.0 

33.01 

7A4-6 

3 

1000 

402 

3105 

2195 

5221 

5711 

25.2 

22.10 


4 

402 

1166 

2195 

2130 

5711 

5630 

27.1 

29.25 

7A4-7 

2 

292 

990 

1711 

3425 

5016 

6446 

9.5 

10.78 


3 

990 

316 

3425 

1838 

6446 

5800 

26.1 

21.68 


4 

316 

1010 

1838 

2120 

5800 

6427 

26.1 

21.90 


t)iff. 


Sec. 

+0 . 3 
-0.7 
+ 0.1 
+1.4 
-3.0 
-4.7 
-1.4 
- 1.6 
-0.7 
-0.7 
+0.9 
+0.4 
-1.3 
-3.1 
- 1.2 

+ 6.0 
-7.7 
-4.2 
+0.4 
+ 3.3 
-3.8 
- 2.1 
+ 1.0 
+ 2.8 
-4.3 
- 2.0 
+ 3.1 
- 2.2 
-1.3 
+4.7 
+4.2 
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TABLE 5 . 

Comparison for Series 407 & 47A 


128 

328 

138 

519 

151 

304 

169 

489 

176 

334 

190 

579 

293 

493 

290 

122 

380 

133 

511 

675 

433 

160 

354 

111 

506 

220 

530 

205 

1034 

213 

721 

1843 

343 

733 

353 

506 

1029 

582 


169 

489 

853 

334 

190 

579 

916 

493 

290 

1004 

380 

133 

511 

1316 

1221 

1514 

354 

111 

506 

1022 

530 

205 

1034 

213 

721 

1843 

712 

733 

353 

1094 

1029 

582 

1518 


1672 

2643 

226C 

3420 

2643 

2459 

3420 

2507 

2459 

2248 

2507 

3135 

2248 

2482 

3135 

3522 

1761 

2845 

3372 

4490 

2845 

2983 

4490 

3220 

2983 

2449 

3220 

4272 

2449 

2618 

4272 

4441 

2207 

2707 

5306 

4887 

2707 

2553 

4887 

4628 

2553 

2176 

4628 

4498 

2176 

2399 

4498 

4515 

2551 

2956 

4262 

5665 

2956 

2887 

5665 

6534 

2887 

2667 

6534 

6222 

1093 

2875 

6202 

5921 

2875 

2592 

5921 

6670 

2592 

2340 

6670 

4515 

2340 

2641 

4515 

5963 

2120 

2573 

3126 

5664 

874 

2601 

4948 

6746 

1441 

2278 

1385 

1988 

2278 

1959 

1988 

1992 

1959 

2079 

1992 

2019 

2079 

2398 

2019 

2302 

1839 

2342 

2384 

3490 

2342 

2316 

3490 

3518 

2316 

3087 

3518 

5897 

3087 

2645 

5897 

4595 

2645 

2439 

4595 

4995 

2439 

2912 

4995 

7472 

2912 

2498 

7472 

6376 

* 

2382 

5698 

7249 

2382 

2041 

7249 

8339 

2041 

1783 

8339 

6614 

2364 

2782 

6388 

8506 

2782 

2586 

8506 

9153 

2586 

2578 

9153 

8381 



22 

.3 

46 

.5 

42 

.6 

30 

.0 

29 

.1 

52 

.3 

54 

.3 

34 

.8 

20 

.0 

29 

.7 

40, 

.1 

28, 

.1 

46, 

.7 

48. 

.9 

32. 

.0 

17. 

.6 

13. 

,6 


41.6 

33.0 

60.0 
32.1 

19.4 
31.9 

39.0 

27.1 

48.3 

28.1 

14.5 

29.1 

39.7 

43.6 

18.4 

24.2 

38.8 

22.3 


43.14 

56.90 

31.63 

22.18 

48.22 

48.78 

31.29 

31.01 

50.50 
55.04 

32.51 
36.03 
33.75 
51.18 
26.11 
50.92 

64.67 

39.67 
17.54 

7.34 

47.20 


36.07 -3.1 
72.61 -12.6 
35.89 -3.8 
17.40 +2.0 
30.13 +1.8 
55.84-16.8 
23.78 +3.7 

51.82 -3.5 
28.67 -0.6 
13.80 +0.7 
31.09 -2.0 

47.86 
18.36 
26.31 

38.83 + 
23.72 1 












TABLE 6 . 

Comparison for Series 




144 

756 

141 

291 

252 

1064 

311 

665 

341 

872 

372 

812 

199 

987 

461 

283 


756 

141 

291 

1274 

1064 

311 

665 

1846 

872 

372 

812 

1844 

987 

274 

283 

663 


1014 

3658 

1458 

1762 

1419 

2943 

2230 

2056 

1682 

2434 

1937 

1895 

2128 

2763 

2551 

2308 


3658 

1458 

1762 

3685 

2943 

2230 

2056 

3197 

2434 

1937 

1895 

2558 

2763 

2401 

2308 

2380 


3482 

4337 

3778 

853 

1623 

796 

6282 

6290 

8437 

5667 

9604 

7395 

15008 

8304 

6176 

9252 

1407 

2688 

933 

2184 

2390 

3487 

2668 

2913 

3618 

4240 

3183 

3461 

3361 

4759 

4316 

3771 


4337 

3778 

5616 

1623 

796 

1663 

6290 

8437 

5667 

5631 

7395 

15008 

6537 

6176 

9252 

8190 

2688 

933 

2184 

2778 

3487 

2668 

2913 

3883 

4240 

3183 

3461 

4105 

4759 

3087 

3771 

3599 


WBl 


iff. 

Exp. I The. 1 


Sec. I Sec. Sec. 


169.4172.41 -3.0 

243.0 3 L 7. 55 -74.6 
93.1100.19 -7.1 


116.1 121.96 
174.3 151.31 
208.0 218.57 
71.6 70.63 


4.70 

6.92 

3.04 

2.74 

3.34 

8.58 

3.32 

1.82 

4.40 

6.82 

2.76 

2.08 

2.82 

8.22 

9.88 

4.84 


-5.9 



+0.7 

+ 1.1 

- 0.6 



No value 
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source was much more stable since it did not have to under- 
go such large changes in heat flux. On a percentage basis 
the longer times were found to compare better. It is point- 
ed out that even with the 0.800 criterion some of the times 
are still very small. It is also noted that the differences 
for these smaller times (904 series) are mostly positive. 
Since the above-mentioned early time accuracy problems would 
tend to cause errors in this direction the effect may still 
have been present. In Table 6 it can be seen that the com- 
parisons are made for only seven of the 405 series experi- 
ments. For the other tests shown there the end conductances 
could not be calculated accurately because the temperature 
drops were too small. This is unfortunate, but nothing can 
be done about it. The limitations of the experimental accu- 
racy (discussed in Section V) prevent any reasonable esti- 
mate from being made. As an example, for a calculated AT 
of 2 °f , a possible error of ± 1°F makes it impossible to 
evaluate the contact conductance within 100%. For those 
runs in which this occurred no attempt was made to guess at 
the conductance values. Although an estimate was made for 
the hot end conductance on two of the phase 1 runs of the 
405 series it could be seen from the theoretical work that 
the effect was small. However, in the present cases the 



182 


theoretical solution shows that the time to approach steady 
state depends on the difference between the initial and fi- 
nal values at each contact. The possible consequences of 
estimating both or even one conductance value are obvious. 
Thus only those runs for which calculated values were avail- 
able were used in the comparison. 

It is difficult to evaluate the comparisons in the 
tabular form. For this reason the data were plotted in the 
form shown in Figure 35. This figure is a plot of the ex- 
perimental values of the time to approach steady state ver- 
sus the corresponding theoretical values. Each series is 
represented by a different plotting symbol as shown on the 
plot. For a perfect correlation each point would lie on a 
45 line passing through the points (1,1), (10,10), etc. As 

a means of evaluating the results two bands are drawn on the 
plot representing an agreement of ± 10% and ± 20 % based on 
the theoretical. There are a total of 92 points plotted in 
Figure 35. Of these only 18 lie outside the 20% band, and 
8 of that 18 are on the 904 series, for which there is some 
doubt about the accuracy. Of the 92 total points, 48 of 
them lie within the 10% band. Although this comparison is 
indirect to the intended purposes, i.e., it could not be 
made in the form of a plot as a function of interface con- 


1 


10 

Theoretical Time - seconds 


100 


Fig. 35. — Time to Reach Steady State 
for All Series-Phases 2, 3, 4, 5 and 6. 
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tact conductance, it is still believed to be a representa- 
tive one. 

In view of the above comparison it is the opinion of 
the writer that the correlation between the experimental 
behavior and the simple, constant property theoretical pre- 
dictions for this second type of transient is very good. 

It is believed that these data are the first data on the 
transient thermal response of composite metals to be re- 
ported. 


VII. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 

Solutions have been derived for a class of boundary 
value problems for the time -dependent temperature distribu- 
tion in a two layer, composite slab with contact resistance 
at the interface and contact or convective resistance on the 
outer boundaries. These solutions represent an original 
contribution to the field of conduction heat transfer. 

The results of a limited parametric computer study 
using these solutions has been presented. This study in- 
cludes a set of dimensionless correlations of an arbitrarily 
defined time to approach steady state and a discussion of 
some transient thermal characteristics which were observed 
for these types of systems. 

To evaluate the usefulness of the above solutions for 
predicting transient response of real systems an experimen- 
tal program was carried out. The experimental work consist- 
ed of measuring the transient temperature distributions in a 
series of test samples when subjected to thermal transients 
which approximated the boundary conditions used in the theo- 
retical solutions. Test samples of aluminum, Armco-Iron, 


185 


186 


and stainless steel were used to provide a large range of 
thermal properties. The test data provided by this experi- 
mental program were compared to the theoretical predictions. 
These data also represent an original contribution to con- 
duction heat transfer. 

For the case of a uniform initial temperature distri- 
bution and a sudden increase in one boundary medium temper- 
ature, the results were compared on the basis of time to ap- 
proach steady state and the time of occurrence of maximum 
overshoot as functions of interface contact conductance. In 
general the agreement was found to be good. For those cases 
in which a noted discrepancey in the experimental boundary 
condition was believed to have only a small effect the agree- 
ment was very good — at least within 10%. 

The results of a second set of experiments for the 
case of an initial steady state heat flow condition followed 
by a sudden change in contact conductance were also compared 
to the corresponding theoretical case. These results could 
only be compared on a basis of time to reach steady state 
for each individual test made because of the large number of 
variables involved. These comparisons showed that less than 
20% of the experimental data differed from the theoretical 
by more than 20%, and that over 50% of the data was within 
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10% of the theoretical. 

On the basis of the above comparisons it is concluded 
that the theoretical solutions presented could be used to 
predict the transient response of systems to which they are 
3-PPli ca t>l e to an accuracy sufficient for most engineering 
purposes. That is, with these solutions properly program- 
med on a digital computer, a knowledge of the required ther- 
mal properties and an estimate of the contact conductance 
from the existing literature, a design engineer could pre- 
dict the transient response of a given system. Thus used an 
accuracy at least as good as his knowledge of the contact 
conductance could be expected. 

Recommendations For Future Work 

recommended that further experimental work be 
performed to provide additional verification of the above 
conclusions. It is believed that emphasis should be placed 
on the phase 2, 3, etc. experiments since these type of 
changes are representative of what happens when a space ve- 
hicle enters a planetary atmosphere. The following recom- 
mendations regarding future experimental work all involve 
improvement in procedures and technique to achieve better 
accuracy. 
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1. Steps should be taken to improve the stability of 
the source block temperature. While a larger 
source block with a larger condensing cavity 
should provide some reduction of this problem, 
consideration should be given to a fast response, 
thermostatically controlled electrical heat source. 

2. For greater accuracy in the temperature profiles a 
larger number of thermocouples should be used, and 
the thermocouples near the boundaries should be 
placed as close to the boundary planes as pos- 
sible. In order to achieve the latter it will be 
necessary to find a better method of contact sur- 
face preparation. The distortions of temperature 
P r °files near the contact boundaries noted in the 
present study and by others [58,92] are the re- 
sult of flatness deviations of the contact sur- 
faces. Great care should be exercised in the sur- 
face preparation to achieve a high degree of flat- 
ness. Since rougher surfaces would also help this 
problem it is recommended that the surfaces be 
ground to optical flatness and then blasted to ob- 
tain non-oriented roughness as discussed by Henry 
[114] . 
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3. It is also recommended that the thermocouple cal- 
ibration procedure be changed to reduce the amount 
of time required for calibration. With the test 
specimens in place in the test fixture, the tem- 
perature distributions could be varied over the 
range of interest by varying the contact pressure 
and/or the source temperature. By use of a multi- 
connected switch all thermocouples could be cali- 
brated with a single potentiometer in a manner 
similar to the method described in section V. 

This method would achieve the same accuracy for 
temperature measurements, but it has the advan- 
tages of being simpler and faster, and would con- 
centrate the calibration points of each thermocou- 
ple into the range over which it would be used 
during the experiments. 

4. Experimental work on thinner specimens would be 
desirable, and consideration should be given to 
this possibility. 

Since the results of the present study are encouraging 
it would seem desirable to expend further effort in the the- 
oretical work. It is recommended that this work be directed 
toward attempting to establish general correlations of the 
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temperature time distributions in a form similar to the 
well-known work of Schneider [251] . 
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P c -contact pressure $ ss =time to reach steady state 

h c -interface conductance 9 mo =time of maximum overshoot 

h l~ hot end " Primes indicate adjusted times 

h 2 = cold end " for h,=h«-h * „ i 

^ ror n^ n 2 "ad. *=No value. 
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APPENDIX B-2 


PHASES 2,3, 4 AND 5 RESULTS 





CONTACT CONDUCTANCES 


Contact 





(Btu/hr-ft 2 - 

°F) 


Pressure 

9 CC 




Hot 

End 

Cold 

End 

In' 1 


ss 

Run 

Ph. 

In' 1 

Fin. 

In ' 1 

Fin. 

In' 1 

Fin. 

psi. 

IS9 

sec. 

304-5 

2 

638 

2319 


4900 

3402 

4500 

35 

EE 

■B 


3 

2319 

1238 


4169 

4500 

4195 

171 

wH 

KH 


4 

1238 

1975 

4169 

4366 

4195 

4305 

35 

■M 

mm 

304-6 

2 

199 

1707 

1200 

4444 

1838 

4366 

11 

105 

9.5 


3 

1707 

446 

4444 

1993 

4366 

2923 

105 

11 

10.7 


4 

446 

1092 

1993 

2345 

2923 

3738 

11 

11 

11.1 


5 

1092 

2441 

2345 

4220 

3738 

4606 

11 

104 

9.2 

304-7 

2 

389 

1893 

1972 

4932 

3225 

4747 

21 

132 

nn 


3 

1893 

907 

4932 

2818 

4747 

3825 

132 

21 

eh ■ 


4 

907 

1549 

2818 

3510 

3825 

4077 

21 

21 

9.7 


5 

1549 

2637 

3510 

4745 

4077 

4917 

21 

132 

8.7 

304-8 

2 

1111 

2987 

4571 

5925 

4904 

5679 

50 

214 

7.2 


3 

2987 

1657 

5925 

4697 

5679 

5400 

214 

50 

9.2 


4 

1657 

2505 

4697 

5063 

5400 

5431 

50 

50 

6.1 


5 

2505 

3553 

5063 

5776 

5431 

6007 

33 

164 

4.4 

7A4-3 

2 

102 

640 

952 

4870 

1063 

3408 

22 

140 

27,1 


3 

640 

121 

4870 

2123 

3408 

818 

140 

23 

53.0 


4 

121 

408 

2123 

2399 

818 

1077 

23 

22 

32.9 


5 

408 

1351 

2399 

4022 

1077 

3716 

22 

140 

37.6 

7A4-4 

2 

117 

868 

1184 

3357 

1127 

4129 

36 

171 

28.1 


3 

868 

146 

3357 

1771 

4129 

993 

171 

36 

53.2 


4 

146 

498 

1771 

1747 

993 

1297 

36 

36 

35.8 


5 

498 

2032 

1747 

3763 

1297 

4505 

36 

205 

24.2 

7A4-5 

2 

196 

686 

1457 

2417 

2998 

5101 

55 

141 

25.1 


3 

686 

234 

2417 

1752 

5101 

3481 

141 

55 

38.8 


4 

234 

808 

1752 

1725 

3481 

4055 

55 

55 

31.0 

7A4-6 

3 

1000 

402 

3105 

2195 

5221 

5711 

193 

81 

25.2 


4 

402 

1166 

2195 

2130 

5711 

5630 

81 

81 

27.1 

7A4-7 

2 

292 

990 

1711 

3425 

5016 

6446 

71 

206 

9.5 


3 

990 

316 

3425 

1838 

6446 

5800 

206 

71 

26.1 


4 

316 

1010 

1838 

2120 

5800 

6427 

71 

71 

26.1 


&ss = Time to reach steady state. 
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PHASES 

APPENDIX B-3 
2, 3, 4, 5 AND 6 RESULTS 






Contact Conductances 

Contact 





(Btu/hr 

-ft 2 - 

■°F) 


Pressure 




1 Interface 

Hot 

End 

Cold 

End 

BIB 

IVTflV 

CT SS 

Run 

Ph. 

In ’ 1 

Fin. 

In ' 1 

Fin. 

In ' 1 

I13B1 

IffiB 

IIHP 

msm 

407-1 

2 

128 

328 

1672 

2643 

2260 

3420 

K 

IK 

40.7 


3 

328 

138 

2643 

2459 

2459 

2507 

HU 

32 

46.5 


4 

138 

519 

2459 

2248 

2507 

3135 

32 

32 

27.2 


5 

519 

1003 

2248 

2482 

3135 

3522 

32 

89 

22.3 

407-2 

2 

151 

304 

1761 

2845 

3372 

4490 

47 

118 

46.5 


3 

304 

169 

2845 

2983 

4490 

3220 

118 

47 

42.6 


4 

169 

489 

2983 

2449 

3220 

4272 

47 

47 

30.0 


5 

489 

853 

2449 

2618 

4272 

4441 

47 

118 

29.1 

407-3 

2 

176 

334 

2207 

2707 

5306 

4887 

64 

150 

52.3 


3 

334 

190 

2707 

2553 

4887 

4628 

150 

63 

54.3 


4 

190 

579 

2553 

2176 

4628 

4498 

63 

63 

34.8 


5 

579 

916 

2176 

2399 

4498 

4515 

63 

149 

20.0 

407-4 

2 

293 

493 

2551 

2956 

4262 

5665 

78 

180 

29.7 


3 

493 

290 

2956 

2887 

5665 

6534 

180 

78 

40.1 


4 

290 

1004 

2887 

2667 

6534 

6222 

78 

78 

28.1 

407-5 

2 

122 

380 

1093 

2875 

6202 

5921 

22 

72 

46.7 


3 

380 

133 

2875 

2592 

5921 

6670 

72 

22 

48.9 


4 

133 

511 

2592 

2340 

6670 

4515 

22 

22 

32.0 


5 

511 

1316 

2340 

2641 

4515 

5963 

22 

72 

17.6 

407-6 

2 

675 

1221 

2120 

2573 

3126 

5664 

29 

86 

13.6 

407-7 

2 

433 

1514 

874 

2601 

4948 

6746 

5 

86 

41.6 

47A-1 

2 

160 

354 

1441 

2278 

1385 

1988 

14 

58 

33.0 


3 

354 

111 

2278 

1959 

1988 

1992 

58 

14 

60.0 


4 

111 

506 

19 59 

2079 

1992 

2019 

14 

14 

32.1 


5 

506 

1022 

2079 

2398 

2019 

2302 

14 

57 

19.4 

47A-2 

2 

220 

530 

1839 

2342 

2384 

3490 

30 

85 

31.9 


3 

530 

205 

2342 

2316 

3490 

3518 

85 

30 

39.0 


2A 

205 

1034 

2316 

3087 

3518 

5897 

30 

173 

27.1 


3A 

1034 

213 

3087 

2645 

5897 

4595 

173 

30 

48.3 


4 

213 

721 

2645 

2439 

4595 

4995 

30 

30 

28.1 


5 

721 

1843 

2439 

2912 

4995 

7472 

30 

172 

14.5 


6 

1843 

712 

2912 

2498 

7472 

6376 

172 

30 

29.1 

47A-3 

2 

343 

733 

* 

2382 

5698 

7249 

53 

131 

39.7 


3 

733 

353 

2382 

2041 

7249 

8339 

131 

53 

43.6 


4 

353 

1094 

2041 

1783 

8339 

6614 

53 

53 

18.4 

47A-4 

2 

506 

1029 

2364 

2782 

6388 

8506 

76 

178 

24.2 


3 

1029 

582 

2782 

2586 

8506 

9153 

178 

76 

38.8 

A 

4 

582 

1518 

2586 

2578 

9153 

8381 

76 

76 

22.3 


®ss - Time to reach steady state. * = No value. 















APPENDIX B-4 
PHASES 2.3,4 AND 5 RESULTS 


Contact Conductances Contact 

L Btu/hr -ft^ — °F ) Pressure 



144 

756 

141 

291 

252 

1064 

311 

665 

341 

872 

372 

812 

199 

987 

461 

283 


243 

84 

801 

127 

31 

343 

168 

49 

711 

1036 

202 I 

89 

1015 

237 

93 

945 

756 

141 

291 

1274 

1064 

311 

665 

1846 

872 

372 

812 

1844 

987 

274 

283 

663 


* 6744 

848 4766 
4766 1481 
1481 1079 


* 7320 

7 320 10001 
2431 
2026 
1885 


7192 

3658 

1458 

1762 

3685 

2943 

2230 

2056 

3197 


1937 

1895 

2128 

2763 

2551 

2308 


3482 

4337 

3778 

853 

1623 

796 

6282 

6290 

8437 

5667 

9604 

7395 

15008 

8304 

6176 

9252 

1407 

2688 

933 

2184 

2390 

3487 

2668 

2913 

3618 

4240 

3183 

3461 

3361 

4759 

4316 

3771 


4337 

3778 

5616 

1623 

796 

1663 

6290 

8437 

5667 

5631 

7395 

15008 

6537 

6176 

9252 

8190 

2688 

933 

2184 

2778 

3487 

2668 

2913 

3883 

4240 

3183 

3461 

4105 

4759 

3087 

3771 

3599 


30 30 


®ss - Time to reach steady state. 


* = No value. 
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APPENDIX C-l: Main Plate and Rod Details 
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Appendix D-4 

Values of Circuit Elements in the Force Control 
System Electrical Circuit Diagrams (pp. 201-3) 


Resistors 


No. 

Value 

No. 

Value 

No. 

Value 

1 

560 

23 

10K 

45 

47K 

2 

15K 

24 

500 

46 

22K 

3 

2.2K 

25 

560 

47 

470 

4 

3.3K 

26 

560 

48 

560 

5 

IK 

27 

10K 

49 

10K 

6 

220K 

28 

120 

50 

10K 

7 

4.7K 

29 

120 

51 

100 

8 

IK 

30 

125 

52 

510 

9 

10K 

31 

125 

53 

1000 

10 

47K 

32 

2.2K 

54 

100 

11 

470 

33 

500 

55 

56 

12 

100 

34 

2.2K 

56 

500 

13 

2.2K 

35 

4.7K 

57 

560 

14 

2.2K 

36 

IK 

58 

220 

15 

4.7K 

37 

68K 

59 

560 

16 

67K 

38 

10K 

60 

IK 

17 

67K 

39 

220 

61 

2.2K 

18 

4.7K 

40 

IK 

62 

500 

19 

2.2K 

41 

4.7K 

63 

250 

20 

2.2K 

42 

470K 

64 

550 

21 

2 . 2K 

43 

470 

65 

2.5 

22 

2.2K 

44 

100 





Capacitors 



No. 

Value 

No. 

Value 

No. 

Value 

1 

5 

10 

.01 

19 

2500 

2 

10 

11 

5 

20 

10 

3 

.01 

12 

10 

21 

10 

4 

5 

13 

100 

22 

5 

5 

100 

14 

250 

23 

5 

6 

5 

15 

10 

24 

.082 

7 

250 

16 

5 

25 

.082 

8 

5 

17 

.05 



9 

.01 

18 

10 




All values in ohms and microfarads, K=1000. 

AH transformers are Thordason— Meissner model TR— 222. 
All transistors are 2N2923. 

All diodes are 1N626. 
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NOTE ON BIBLIOGRAPHY 

The standard abbreviations have been used for all 
journal articles cited. Following is a list of abbrevia— 
tions used for some reports and translations. 

Reports 


USAEC - United States Atomic Energy Commission. 

UKAEA - United Kingdom Atomic Energy Authority. 

Translations 

RSIC - Redstone Scientific Information Center, Hunts- 
ville, Alabama. 

DSIR - National Lending Library for Science and Tech- 
nology, Boston Spa, Yorkshire, England. 

SLA - Special Libraries Association: Translation 

Center, The John Crerar Library, 35 West 33rd St., 
Chicago, Illinois. 

ASLIB - Association of Special Libraries and Information 

Bureau, 3 Belgrave Square, S.W. 1, London, England. 
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